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SYNOPSIS 


Flow in alluvial channals is often accompanied by 
inf luent/eff luent seepage through the channel bed due to different 
free surface levels in channel and ambient ground water table. 
The problem of quantifying the effect of cross-bed seepage on open 
channel (OC) flow superficially resembles that of distributed 
suction/injection on a turbulent boundary layer (TBL), but the 
usual TBL assumption of constant Vfy) (= injection rate V^) is 
invalid in OC flow. The one-dimensional analysis commonly used in 
OC hydraulics is also unsatisfactory due to non-description of the 
mean flow and stress fields and assumptions on friction slope, S£ . 
In this work streamwise velocity, U, bed shear velocity, U* , and 
energy gradient velocity, a^HS (a^ is gravitational 

acceleration, H is flow depth and S is free surface slope) are 
assumed to be piecewise linear in X, 

U(x,y) = FCy) + X GCy), x = X/Hq, Y = Y/K^; ] 

U*(x) = U* Cl + (3x); a^HS = (1+ ax) j (1) 

It is shown that when S ~ (bed slope boundary layer type 
equations are obtained if IVy/U^I << 1, which yield two coupled 
ODEs involving Ffy) and G(y) . Approximate solution of the ODEs 
for |Vyl<< 1, where v^ = Vy/'^agH^S^, shows that the flow is 
nearly self-similar with G(y) oC FCy) and 

FCy) = in [y*' - y^*^ - 6n^ CVv - t/yo)] C2) 

where n varies directly with v^, y^ is determined using Rubesin's 
criterion of invariant Reynolds number at sublayer-turbulent 
interface, and a~ P — GCy). which resemble 

power laws, approach the log-law identically as v^ — > 0. 
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A aet of flume experimenta with immobile bed waa 
conducted to verify the above model. Hydraulically rough, 
aubcritial OC flow waa subjected to inf luent/ef fluent seepage 
through permeable sand bed. Sand bed thickness was varied with X 
to ensure uniform . Discharge, flow depth and water surface 
profiles were measured. It was found that, within limits of 
experimental error, both the water surface curvature and the 
change in from uniform flow depth (for same discharge) compare 
favourably with analytical results. 

The following results follow from the analysis : 

2 

Ci) U* /CagH^SQ), or Sf/S^ where Sf is friction slope, 

increases/decreases with suet ion/ in j ection in agreement with 
one-dimensional OC analysis and TBL data. 

(ii) a and 3 are comparable, but deviate as |v^| increases. 

(iii) Dimensionless velocity u(y) increases and has a fuller 
profile for suction; conversely for injection. 

(iv) ^(y) sharply increases/decreases near the bed with 

suction/inj ection. 

(v) The above changes magnify with increase in Hq SgH and 

with decrease in relative roughness, Kq/Hq. 

Cvi) Uater surface has concave/ convex curvature for 

suet ion/ in j ection . Flow depth varies as 

H - Ho [1 -a SoX^ {1 + (a / 3 )x)] (3) 

Cvii) The analysis suggests that in duct flow dP/dX 

decreases/increases with suction/injection for same 
discharge . 

(viii) From a physical interpretation of parabolic eddy viscosity 
distribution it is reasoned that the flow is nearly self- 
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preserving in all respects — mean flow, turbulence field 

2 

and coherent structures — if |Vy U/U*l << 1 , and an 

estimate of inner boundary layer length with change in wall 
inputs is suggested. 

(ix) It is suggested that effluent seepage increases flow depth 
and enhances chances of flood waves overflowing stream 
banks, while infiltration may be the cause exponentially 
decreasing bed slope in graded streams. 

Seepage effects were included in an investigation of 
longitudinal ribbons (or stripes/lineations) of lateral grain 
sorting and relief on non-uniform sand beds. Subcritical flow 
experiments with equilibrium sediment transport rates (Qg) were 
conducted with beds of three different sand mixtures SI, S2 and S3 
having median diameters, DgQ = 0.28, 0.52 and 1.14 mm and Kramer's 
uniformity coeffcient of 0.48, 0.35 and 0.18 respectively. 

Flume runs with and without seepage produced ripple 
bedforms on SI sands even at low Q^. For S2 and S3 sands stable, 
uniform streamwise stripes developed at low Qg, which broke down 
into primary ripples/dunes at higher sediment load. The lateral 
stripe spacing, L^, was typically about 1.94H to 2.15H, with 
relief, Ly, of about 1.5 to 3.0 D 50 . and with considerable grain 
sorting (fine grains on ridges and coarse grains in furrows). Qg 
was found to increase/decrease with ef fluent / inf luent seepage but 
was in reasonable agreement with the Einstein-Brown relation. 
Lateral grain sorting and Ly were weakly dependent on Qg ( and 

hence on V^) . The lineations covered a significant region of the 
lower and transition bedform regime as delineated by Simons and 
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Richardson. It was concluded that the appearance of three— 
dinensional 'bedforms at low to moderate Qg is inhibited for very 
coarse sand-fine sand or sand-gravel mixtures, which complements a 
similar finding by Chiew for supercritical flow conditions. 

The flow dynamics governing lineations has been related to 
secondary flow of Prandtl’s second type by several researchers, 
whose modelling of the flow and boundary stresses involve 
assumptions on U(y) and U(y), or on V'^CViZ) and l/^CyiZ)i besides 
other approximations. It is suggested here that an eddy viscosity 
(or mixing length) formulation with other appropriate assumptions 
may be used in. the X-momentum equation to obtain the approximate 
flow field and boundary stresses. 

The major conclusions of the thesis are : 

(i) OC flows over long stretches with low seepage rates are 
almost self-similar. The present results for U(y), T (y) , 
H(x) and S£ can be directly used for known seepage rates in 
unlined canals, alluvial streams and other OC applications 
with clear water flow. 

(ii) Nonuniform or sand-gravel beds develop stripes of lateral 
grain sorting and remain nominally plane at low to moderate 
sediment load in subcritical flow. Lateral sorting and 
secondary flow cells are dynamically related through the 
redistributed stress field. 

(iii) Mixing length or eddy viscosity approach may be gainfully 
applied to OC (and other wall turbulent) flows when 
secondary velocities (V and/or U) are small in comparison to 
the streamwise velocity. 
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CHAPTER 1 


INTRODUCTION 


1.1 General 

Surface water noveinent throufih natural and artificial open 
channels plays a fundamental role in the distribution and 
availability of water for human needs and concerns, such as 
domestic, agricultural and industrial use, aquatic life, sediment 
transportation, navigation, power generation, waste disposal and 
ecology. While the complexity of these processes in the natural 
environment is daunting, considerable understanding of alluvial 
channel mechanics hais been achieved by applying the basic 
principles of fluid notion. Application of these principles to 
open channel flow is mainly governed by the nature of flow 
laminar or turbulent, the surface and body forces acting in a 
given fluid domain, and the characteristics of the free surface 
and the solid boundary. 

Flow in alluvial channels is usually turbulent, the 
instantaneous velocities of small fluid "elements” being random 
and three dimensional even when the mean flow is steady and 
rectilinear. Moreover, in turbulent flow, fluid stresses act 
predominantly through correlated velocity fluctuations rather than 
mean velocity gradients. As velocity fluctuations tend to die out 
in a region very close to the boundary — the viscous sublayer, 
the impulse of these fluctuations contribute to much shear stress 
and pressure fluctuations at the boundary. The process is more 

complex when the boundary is hydraulically rough which is common 

in gravel- and sand-bed channels. In such cases random 

protrusions on the boundary extend beyond the sublayer, 
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effectively destroyinfi it, and affect the near— bed turbulence 
fleneration. The effects of a rou£h boundary, transmitted through 
turbulent eddies, may permeate the entire boundary layer. 

Since fluid particles immediately adjacent to a solid 
boundary tend to adhere to it in both laminar and turbulent flows, 
the flow is bounded by the no-alip and no-penetration conditions 
at a rifiid, impervious surface. However, a boundary composed of 
layers of loose sediment, as often encountered in sedimentary 
hydraulics, is often permeable, which fienerally implies non-zero 
mean and turbulent velocities at the boundary. Even if the mean 
velocity vector at the boundary be known, the corresponding 
fluctuating quantities are difficult to determine. 

Alluvial channel mechanics also involve processes other 
than these encountered in pure fluid mechanics, chiefly because of 
sediment transported by the flow. Sediment is mainly transported 
in two ways — as suspended sediment and as bedload. Suspended 
sediment grains are believed to be sufficiently small as to 
respond directly and continually to small scale velocity 
fluctuations. Hence their effects can be accommodated in the 
equations of motion for turbulent flow. The pre-eminent view 
considers suspended sediment as passive elements or contaminants 
which, while responding to the fluid forces exerted on them, can 
significantly alter flow properties through some form of passive 
resistance. In contrast to suspended sediment, bedload motion 
consists of intermittent sliding, rolling and saltating movements 
of individual grains, which, being discrete events in space and 
time, are not amenable to continuuie analysis. However, by looking 
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at bedload motion as essentially a random, surface phenomenon, it 
is possible to assimilate its effect in bulk or averafie terms. 

Apart from its direct effect on flow properties, sediment 
motion can have an indirect but important effect on the flow. 
This is because sediment motion often leads to the development of 
bedforms. Bedforms are relatively immobile or slow-movin£ sundry 
deposits of sediment on the bed with fairly regular dimensions and 
patterns. Bedforms may be sufficiently large as to decelerate the 
flow and offer considerable form drag in addition to the native 
skin friction of the bed. Their effects on the mean and turbulent 
velocity and stress fields are also significant. 

1.2 Statement of the problem 

The flow in an alluvial channel is often subjected to seepage 
or infiltration through the porous channel boundary. Seepage 
occurs due to piezometric gradient between the channel and the 
ambient water table. Higher ground water table results in 
discharge into the channel (effluent seepage or injection) while a 
lower ground water table causes abstraction from the channel 
(influent seepage or suction). Available information [eg. Uatters 
and Rao, 1971; Brink and Oldenziel, 1974] indicate that seepage 
affects the hydrodynamic forces on the bed and, possibly, the flow 
and stress fields. There may thus be significant changes in bulk 
flow properties such as flow depth, bed shear stress and the 
balance between the flow, sediment discharge and bedforms. 
Correct prediction o£ the channel regixfte for different eeepage 
conditions will necessitate quantitative modellinft of 
suction/injection on loose boundary open channel flow. 
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In addition to seepaae, the flow in alluvial channels is 
commonly associated with bedforms. Uhile most bedforms exert 
considerable pressure dra« on the flow, sand ribbons (stripes or 
lineations) which are flow-aligned structures on the bed, offer 
little additional flow resistance. Consequently, the discharge 
carrying capacity of channels is not much affected by these 
structures in contrast to other types of bedforms such as ripples 
and dunes. Uhile some researches have been conducted in the past 
few decades [eg., Mantz, 1978; McLean, 1981; Allen, 1982; 
Tsujimoto, 1989], the bedform regime (existence criteria) of sand 
ribbons, their mechanics of formation, and the associated 
structuring of the flow remain incompletely understood. 

1.3 Scope of the work 

The present work attempts to study the above problems 
theoretically and experimentally with a view to gain practical 
insight. Earlier studies have shown that secondary velocities 
introduced in the cross-flow plane by seepage and sand ribbons are 
much smaller than the strearowise velocity component almost 
throughout the flow region. The low secondary velocities are 
unlikely to cause significant flow curvature or altered turbulence 
generation processes, and it may be anticipated that no 
fundamentally new assumptions are necessary to analyze the flow. 
The analytical work attempted herein is limited to defining the 
time averaged "velocity and stress fields from the equations of 
motion, and detailed understanding of the turbulent flow structure 
is beyond the scope of this work. The analysis is restricted to 
the case of steady, clear water flow in wide channels. 
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The theoretic&l work has been au£inented by experiments 
designed to verify the analytical results as well as to provide an 
independent data base for dravinfi empirical conclusions. In 
P^-fticular , experiments were needed to quantify the existence 
criteria and physical properties of sedimentary bedforms with 
equilibrium sediment transport. Measurement of the flow structure 
associated with sand ribbons, which has been reported by earlier 
researchers, was not attempted in this work. 

1.4 Thesis presentation 

The thesis has been presented in six chapters. In the 
following chapter the literature on open channel (and other types 
of wall-bounded) flow has been reviewed. The effects of seepage 
on open channel flow and the dynamics of channel flow with sand 
ribbons as reported in the literature have also been reviewed in 
Chapter 2. In Chapter 3 analysis of the mean flow field subjected 
to suction/ in j ection of fluid through bed and in the presence of 
sand ribbons have been separately carried out. The experiments 
that were conducted are described in Chapter 4. In Chapter 5 
experimental results as well as data from literature have been 
analyzed and discussed in association with the analytical models 
of Chapter 3. A brief summary and major conclusions of the thesis 
have been presented in Chapter 6 along with a mention of certain 
issues for potential research ensuing from this study. 



CHAPTER 2 


LITERATURE REVIEW 


In a fluvial channel the hydrodynamic forces acting on the 
boundary, the rates of sediment erosion and deposition, and the 
water and sediment discharges are intricately related. Hence, 
simplifications and idealizations are necessary to study any of 
the component processes. In this chapter literature on the effect 
of suction and injection on channel flow and the mechanics of 
streamwise bedforms has been surveyed. Some of the basic 
assumptions involved in moat studies, which are not specifically 
mentioned in this survey, are the conditions of steady, uniform 
rectilinear flow, steady and uniform sediment transport rate, and 
no long-term channel adjustment or time-dependent environmental 
inputs (such as temperature). The literature has been broadly 
grouped under the following heads: 

* riean, i.e time-average, velocity and boundary forces in open 
channel flow (hereinafter, OC flow). 

* Hydrodynamic effects of seepage inflow/outflow on channel 

bed . ! 

* ITechanics of longitudinal bedforms. 

2.1 Ilean velocity and boundary forces i 

2.1.1 Depth averaged mean velocity ! 

Empirical equations for the depth averaged mean velocity, 
have been widely used in OC hydraulics. Among them the most 

familiar ones are .the Chezy’s and Manning’s equations [Chow, 


1959] : 
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Uav = C : Chezy 2.1 

: Manninc 2.2 

n 

where R is hydraulic radius, S is the energy firadient of flow, and 
C and n and Chezy's and Manning’s respective resistance 
coefficients. For coarse channel material, C and n may be 
estimated from grain size, but in general both coefficients are 
determined empirically to account for a host of channel properties 
[Chow, 1959]. Thus, despite their simple form and ease of use, 
Equations 2.1 and 2.2 are not ideal for making accurate 
predictions or to study the effect of specific variables on the 
flow. A better understanding of the flow process and description 
of mean velocity and stresses have been acheived by drawing upon 
diverse studies on turbulent wall (or wall-bounded) flows since 
the beginning of this century. A brief review of turbulent 
boundary layer, open channel and pipe flows is made here, 
focussing particularly on the mean flow structure. 

2.1.2 Mean velocity in wall flows 

Turbulence in wall flows is characterized by a non-random 
flow field that is apparently related to the local mean velocity 
gradients and stresses, and is governed by the no-slip and no- 
penetration conditions at a rigid, impermeable wall. In terms of 
the mean properties of flow, a turbulent boundary layer 
(herinafter, TBL) consists of the following main divisions 
(Cantwell, 1981; Hinze, 1975]: 

(i) The viscous sublayer is a thin region immediately adjacent to 
the wall and about five wall units, i.e. 5 V /U* , thick, — V 



ia kinematic viscosity of fluid and U* is the wall shear velocity. 
Due to the no-slip condition at the wall the mean velocity in 
sublayer is small as compared to the outer flow, and the shear 
stress is mainly viscous or "laminar” stress. U(y) is well 
described by the laminar flow equation : 

u=y:0<y<5 2.3 

where u = U/U* , y = YU*/?; . U* = / P » ^ being the normal 
distance from the wall, t the wall shear stress and p , the 
density of the fluid. 

(li) The buffer layer, extending from the outer edge of the 
viscous sublayer to the edge of the turbulent region, has 
typically the properties of both adjoining regions, the viscous 
and turbulent stresses being of comparable orders of magnitude in 
this zone. In the buffer layer U(y) is best described as a 
smooth changeover from Equation 2.3 to Equation 2.4 (stated 
below) . 

(ill) The turbulent region covers the boundary layer from 7 ~ 30 
outwards, and predominantly governed by turbulent or Reynolds 
stresses. Significant correlations of velocity fluctuations [see, 
e.g. , Sabot and Comte-Bellot , 1976] clearly indicate a time-mean 
structuring of the flow, which may be described by suitable 
length and time scales of the time-averaged flow (eg. , mixing 
length and wall shear velocity), or by turbulence scales (e.g. 
integral length scale and Taylor microscale), or by scales of 
identifiable coherent structures (eg., the dimensions, lifetimes 
and vorticity of persistent eddies), all of which are related in a 



coinpoai'td flov process. The turbulent region nay be further 
subdivided into two zones, namely: 

a) An inner or lofiarithnic re&ion extending upto about O.IH to 
0.2H, H being the TBL thickness, having the mean velocity 
distribution 
1 

*2 = In y + B 2 4 

X 


where In denotes natural logarithm and X and B are empirical 
constants nearly independent of the flow Reynolds number, Re, 
in canonical TBLs . In the major part of the inner region u 
also varies logarithmically with the outer variable y = Y/H 


rv* 1 ^ 

u = In y + C 

X 


2.5 


C being an empirical constant. Thus the inner region is also 
an overlap region where the velocity scales with both inner and 
outer flow variables. 

(b) An outer turbulent zone, 0.15 < y < 1.0, where u scales on y 
though not logarithmically, the flow turbulence is 
intermittent, and the outer edge undulates along its boundary 
with the irrotational flow region. For constant pressure 
flows the velocity is well described by Hama’s defect law 
[refer Hinze, 1975] : 

Utj - u =9.6 Cl-y)^ 2 . 6 


where is the velocity at the TBL edge. For pressure 

gradient flows ie., 3P/3X^0, where P is the freestream 
pressure and X is the downstream direction, the velocity 


profile is aptly described by Coles’ law of the wake : 


10 


u = --- In y 4 C + U^Cy) 
where = 2 sin^(iry/2) and 

value depends on 3p/ 3X. For 
following empirical relation: 

6 * 9 p n n F 

r = 0.8 (-.j— ---- 4 0.5)°-^^ 


2.7 

r is a wake parameter whose 
P , Uhite [1974] suggested the 


2.8 


w 


where ^ * is the displacement thickness. For nearly 
equilibrium TBL flow in adverse pressure gradient, Schofield 
-,[1981] indicated that the outer region velocity profile can be 
collapsed into a single relation between 0/0^ ^ *• 

the universality of this relation has been questioned by Dengel 
and Fernholz [1990]. 

For turbulent flow in smooth-walled pipes (as also in other 
prismatic ducts), there is neither any intermittency in the flow 
nor any variation of kinematic quantities in the flow direction. 
The sublayer velocity distribution is identical to that in a TBL, 
and the log-law holds in a major part of the turbulent region, 
with some change in the additive constant. Thus in terms of wall 
variables and in defect form. 


^ 1 ^ 

u = In y 4 5 . 5 2.9 

;c 

1 

Uj^ - u = - ---In y 4 0.8 2.10 


where H is the pipe radius. In the core region the velocity 
deviates from the above distributions. Nikuradse had empirically 
shown that in the log-region the velocity can also be approximated 
by a power law, 
u = y" 


2.11 



where Cjj is an empirical constant dependent on n, and 0 < n < 1. 

Typically n decreases with increase in Re so that Equation 2,11 
approaches the log-law asymptotically. 

Equations 2.3 through 2.11 are valid for flow over smooth 
boundaries. For flow over rough surfaces the sublayer flow is 
disrupted while in the turbulent region the velocity is 
depreciated. Based on the early works of Nikuradse [refer 
Schlichting, 1979] and subsequent experiments by many researchers 
[eg., Bayazit, 1976; Stukel, Hopke and Nourmohammadi , 1984] the 

modified log-law over fully rough surfaces 

1 Y - Ay^ 

Uj. = C ) + Bi 2.12 

* k 


where subscript r denotes rough boundary. Kg is Nikuradse 's 
roughness scale, Bj^ is a function of Kg, and A Y^j is a downward 
shift of apparent velocity origin. On comparison with Equation 
2.4, the velocity shift A'u = u^ - u is 

r.» ^ r.j 

A u = Bj^ - B In kg 2.13 


Physically three flow regimes are distinguishable depending on the 
magnitude of "kg , viz. 

* Hydrodynamically smooth regime, Icg < 5, where the sublayer 

** «v> 

flow is virtually undisturbed; and = B = 5.5. 

* Transition regime, 5 < kg < 70 , where the sublayer is 

partially destroyed; Bj^ increases monotonically to about 9.6 
at kg ~ 12 before decreasing to about 8.5 at kg ~ 70. 

* Rough flow regime, kg > 70, where the sublayer is effectively 
destroyed; Bj^ 21 ® • 5 • 
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In contrast to Nikuradse’s three dimensional k-type 
roufihness , Perry, Schofield and Joubert £1969] distinfiuished two- 
dimensional d-type transverse roufihness elements where the 
rounhness scale determining the velocity shift — logarithmic 
asymptote is proportional to H rather than the roughness 

height. Perry et al postulated that the vortices formed behind d- 
type roughness elements are stable and do not shed periodically. 
Their empirical results have been confirmed by later studies, eg. 
Stukel et al [1984]. 

2.1.3 Turbulence structure in wall flows 

Spatial and temporal correlations of different flow 
variables, referred to as coherent motion [Robinson, 1991], 
indicate a stable mechanism in the maintenance of Reynolds 
stresses and the production, convection, diffusion and dissipation 
of turbulence in wall flows. Najor results of coherent flow 
studies have been extensively discussed and reviewed in 
literature, in particular by Cantwell [1981] and Robinson [1991], 
and are discussed here briefly and selectively in order to 
highlight some complexities of near wall flow vis-a-vis the mean 
flow description of the previous section. In terms of persistent 
flow vortices with a directional preference, the outer region of a 
TBL is characterized by weak spanwise vorticity of relatively long 
lifetimes. The bulges and vorticity at the outer edge of TBL 
enable entrainment of freestream fluid into the boundary layer. 

At the other extreme, the quiescent sublayer is densely 
covered with high and low speed streaks ie, alternate elongated 
regions of high- and low-str eamwise velocity, with spanwise 



wavelength of about 100 ^ /U* . The streaky sublayer structure is 
believed to be associated with contra-rotatine pairs of streamwise 
vortices, which are probably the stretched limbs of U- or V-shaped 
horseshoe or hairpin vortices [Blackwelder and Eckelman, 1979; 
Gupta et al . , 1971; Kline et al, 1967]. As compared to outer 
re£ion eddies, the lifetime of sublayer eddies is rather short due 
to quasi-periodic disturbances (the burst cycle) when high speed 
fluid from the turbulent region invade the sublayer, followed by 
lift-up and ejection of low speed streaks into the turbulent 
region. As of now there appears to be no clear consensus about 
whether the period of turbulence bursts scales with inner or outer 
flow variables [Shah and Antonia, 1989]. Moreover, though 
sublayer vortices are believed to break down or to merge with 
neighbouring eddies during bursts [Perry, Henbest and Chong, 1986; 
Smith and Metzler, 1983; Uark and Nagib, 1991], they are also 
reported to be significantly present in most of the TBL [Head and 
Bandyopadhyay , 1981]. 

Near-wall measurements indicate high correlations between 
instantaneous wall pressure, shear stress and velocity components. 
High wall shear is associated with high spanwise vorticity near 
the wall, the rms wall shear fluctuation, ”^v,rms about 0.4 

y [Alfredsson, Johansson, Haritonidis and Eckelman, 1988]. High 
wall pressure, associated with the passage of inclined shear 
layers in the inner turbulent region, have rms values, P'y^rms' 
about 2 to 4 [Hinze, 1975; Uillmarth, 1975]. This value is 
much higher than '’'y^rms Probably due to significantly higher 
contribution of high frequency Py as compared to ^ 

Kuzentsov and Sysoyev, 1983]. 


[ Yef imtsov 



For turbulent flow over roush surfaces, where the sublayer 
and buffer layer may be fully disrupted, the wall-normal and 
spanwise velocity components, V' and U' , are an the whole reduced 
in the loflarithmic zone but do not affect the wake region. 
Apparently, roughness elements increase the mass and momentum 
exchange in the near-wall region, but do not alter the nature of 
flow, the mode of turbulence production and Reynolds stress 
maintenance (Bushnell and McGinley, 1989; Grass, 1971; 
Nourmohammadi , et al , 1985; Sabot et al , 1977]. For densely 

packed roughness elements both ^w increase with 

increase in , but the ratio decreases [Blake, 1970]. 

2.1.4 Naan velocity distribution in OC flow 

The mean velocity profile in turbulent open channel flow 
has been observed to be more or less logarithimic at sufficient 
distance away from the wall. The effect of bed roughness is also 
similar to that for the inner region of TBL flow. However, due to 
natural sediment coming in varying compositions of different 
shapes and sizes, the estimation of poses considerable 

difficulty in natural channels. David [1980], Garde and Raju 
[1985] and Raudkivi [1990] have surveyed various formulae 
recommended for estimating , a few of which may be cited here : 
Ackers and White : Kg = 1.25 D 35 

Einstein and Barbarossa : Kg = 0^5 

Engelund and Hansen : Kg = 2.15 

Hey : Kg = 3 . 5 D 34 

where D is grain diameter, and its suffix denotes the percentage 
of finer grains in the mixture. It is evident that the criteria 
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of Equation 2.15 can give widely divergent results for the sane 
sediment bed. Considering such discrepancy to be due to 
representation of a heterogeneous sediment mixture by just one 
statistical size, Mittal [1977], Mittal [19 8], David [1980] and 
Singh [1981] conducted extensive measurements with fixed beds of 
non-uniform sand and gravel, in which D^q was varied from 0.14 to 
3.0 mm combined with variation of the standard deviation of grain 
sizes, (5jj. Based on their data Gangadharaiah , David and Roy 
[1991] found that K^/D^q can be related quite accurately to 
( “ Their plot shows that Kq/Djq increases roonotonically 

from 1 to about 4 . 8 at Pd ~ 70 before decreasing to about 1.7 at 

rk» 

= 1000- The peaking of the curve at = 70 shows similarity 

of the effect of with that of kg as it coincides with the value 

fO 

of kg when a flow becomes fully rough. A suitable empirical fit 
of Gangadharaiah et al’s plot may be 



r I 

1 + 3.75 exp -1.3 llogiQ 



2.16 


The eatimatlon of bed roughness with sediment transport is 
difficult, but it is known that Kg may change drastically from 
fixed bed values. In particular, Uilson [1989] has shown that 
with bedload motion at high shear stresses K^ is about half the 
shear layer thickness of moving bedload. Uilson also gave the 
following empirical relation 


1th 



ag ( Pa - P ) D 


2.17a 


2.17b 


where ia the aediotent mobility number or Shlelda dimenaionleaa 

ahear atreaa parameter, la acceleration due to gravity, and 
and P are the denaltlea o£ aedlment and water reapectlvely. 

Uhlle moat OC atudlea Involve aubcrltlcal flowa, a few 
atudiea of aupercritlcal flow velocltlea on plane beda have 
indicated aome devatlon from the canonical diatrlbution . In a aet 
of experlmenta involving mean and fluctuating velocity 
meaaurementa Tomlnaga and Nezu [1992] have recently reported that 
though the velocity diatrlbution followa the log-law In the main, 
the additive conatant of the law decreaaea aubatantatlal ly with 
Increaae in Froude number, the decreaae being more algnlfieant for 
rough bed flowa. 

Several OC velocity meaaurementa for aubcrltlcal flowa have 
alao indicated that u de v| atea from the log-law at dlatancea 
beyond about 0.1 5H from the bed, eg. Bronco and Parthenaidea 
[1971], Cardoao, Graf and Guat [1989], Coleman [1981], Finley 
[refer Townaend, 1976], Nezu and Rodl [1986], and Zlppe and Graf 
[1983]. Coleman propoaed to explain the devlatlona by uae of 
Colea' wake law with the modified value of T = 0.191 for clear 
water flow. Noreover, in aome of the data of Cardoao et al and of 
Nezu and Rodl the velocity profile deviatea from the auggeated 
wake law near the free aurface where the flow la apparently 
retarding. Retarding flow near free aurface la. In fact, not 
uncommon In OC flow, particularly when the bed ia hydraulically 
rough [Chow, 1959], and may be attributed to aecondary flow 
[Cardoao et al ] or wind ahear at the free aurface. The latter 
poaalblllty la normally dlacounted in OC flow analyala, but may 
not be entirely without importance. Several atudlea, eg. Raahldl 
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atid Banerjee [1990], Street [1979] and Ueda et al [1977] have 
indicated the poaaibility of free aurface ahear producing aimilar 
effecta aa that of a moving wall, and thia analogy haa been 
recommended for modelling of aurface ahear and near aurface flow 
[ASCE Taah Committee, 1988]. 

In addition to log-law and wake law, power lawa have alao 
been employed for describing OC velocity distributions . In recent 
times Chen [1991] haa demonstrated the validity of power lawa in 
OC flow. He recommended the use of the 1/7-th power law (ie, n = 
1/7 in Equation 2.11) for hydraulically smooth flows, and the 1/4- 
th power law for hydraulically rough flowa. Coleman [1991] has 
alao expressed preference for power lawa over the log-law. 

2.1.S Analytical derivation of mean velocity and bed stresses in 
OC flow 

The two dimensional analysis of flow in wide open channels, 
in its simplest form, uses the continuity and momentum equations, 
and follows the procedure of Reynolds stress modelling in the 
inner region of a TBL. Though several analytical methods exist 
for obtaining the log-law in the inner region, eg. , the overlap 
layer criterion of Millikan [refer Hinze, 1975] and the vorticity 
transport hypothesis of Taylor [refer Schlichting, 1979], 
Prandtl'a phenomenological approach is the moat widely used one 
due to its insight and simplicity of use. According to Prandtl’a 
hypothesis for TBL flow : 

(1) The inner turbulent region is of constant shear stress (='^w) 
(ii) 0’ and V’ are proportional to the local mean vorticity as 
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with - P U’V’ being the time- averaged stress. Uith T = the 
logarithmic velocity distribution is readily obtained. Based on 
Prandtl’s approach, alternative formulations of mixing length or 
eddy viscosity have been proposed by many researchers for other 
types of wall flows and for non-canonical TBL flows [see eg. , 
Galbraith and Head, 1975], 

In order to derive a single velocity distribution for the 
entire inner region including the sublayer and buffer layer, 
Prandtl’s mixing length formulation has been modified or adapted 
by several researchers, eg., by Rotta and by Van Driest 
respectively as [see Hinze, 1975] 


^ , 9 I • 

= V + X ^ I I 

I 3 Y 1 

, , y I 3U I 

and V. = ~ expC- )] 1 I 

A I 3Y I 


2.21 


2.22 


where A is a constant (= 26). More recently Haritonidis [1989] 
has argued that, instead of Equation 2.18, V’ should be determined 
from the frequency of turbulence bursts. Thus he derived 
= \y Y 


2.23 



where is an empirical constant. With certain assumptions and 
qualifications his results show good agreement with the 

fs> 

distribution of u and T in the near wall region of both ordinary 
and drag reducing flows. 

Fully developed open channel flow of high aspect ratio is 
also amenable to the mixing length approach of Prandtl. However 
since there is virtually no constant stress zone in OC flow — 

T increases linearly from the water surface to the bed — 
alternative forms of or Prandtl ’s 1 have been suggested, since 

Prandtl ’s hypothesis is, in any case, intuitive and not based on 
constitutive relations [see eg, Rijn, 1984; Umeyama and Gerritsen, 
1992]. Two such commonly employed formulae cited below are the 
modified mixing length formula and the parabolic eddy viscosity 
formula respectively : 

2 2 ' ' 

\ Y Ci-y) I 1 2.24 

I 9Y I 

and 2.25 

Both equations give the log-law in steady, uniform OC flow for 

(1-y)- 

As mentioned in the previous section the inclusion of a 
wake component in the mean velocity has been suggested in several 
studies. Thus Vedula and Achanta [19SS] included Hama's defect 
law for y > 0.15 in their analysis and found good agreement of 
observed and computed mean velocity distribution and bed shear 
stress for both pipe flow and OC flow data. Willis [1985] also 
Included a wake component to obtain fair agreement with outer 
channel flow data, while for flow in the near-bed region a 
combination of fluid and eddy viscosities yielded good agreement- 



The simple phenomenological approach o£ Prandtl does not 
often predict the velocity distribution accurately in non- 
canonical vail flows. In the past few decades higher order 
turbulence closure models for the Reynolds stresses have been 
increasingly used. These models, which are more complex than the 
one-equation gradient transport formulae, usually need more 
empirical data (for validation/calibration) and computational 
effort. They have been reviewed in fair detail elsewhere [ASCE 
Task Committee, 1988; Speziale, 1991]. 

The time-averaged pressure, P^, and shear stress acting on 
the channel bed is easily computed from the depth-integrated 
momentum equations in X and Y directions [Chow, 1959]. For 
channel with small slope, S, and V’ = 0 at the free surface. 


Pw 

= /H 



2.26 

and Ty 

= y HS 



2.27 

where "Y is 

the unit 

weight of 

water. The mean 

shear stress 

acting an 

individual 

bed grains 

may, however, be 

significantly 


different from that of Equation 2.27 depending on the presence or 
absence of bedforms, and on the arrangement, concentration and 
sorting of grains on the bed surface. Surface grains are also 
subjected to considerable pressure and shear stress fluctuations 
as discussed in Section 2.1.3 and in Vanoni [1977]. 

In mobile bed channels the bed shear stresses and channel 
resistance are not merely dependent on the roughness of bed 
grains. They also depend on the amount and type of sediment load 
(i.e., suspended sediment load or bedload) and other possible 
effects such as sidewalls, bank vegetation and free surface shear 



produced by vinda . Empirical and aemi— empirical formulae have 
been developed by aeveral researchers to account for the effect of 
sediment motion sn OC flow eg, by Einstein and Barbarossa, 
Engeland and Hansen, and Ackers and Uhite [refer Ackers, 1983; 
Vanoni, 1977]. However, it nay be mentioned that due to lack of 
in-depth understanding of the processes , there are wide 
differences in the physical concepts used as well as the sediment 
load predicted by different formulae [ McBean and Al-Nassiri , 
1988; Samaga, Raju and Garde, 1986], 


2.2 Flow over permeable boundary with distributed suction or 
injection 
2.2.1 General 


Turbulent flow over a permeable boundary accompanied by 
suction or blowing of fluid through the boundary may occur in 
diverse natural situations. For example, in rivers and canals 
flowing in granular material, seepage inflow from the channel to 
the bed (influent seepage, infiltration or suction) or seepage 
outflow from the bed into the channel (effluent seepage or 
injection) is a common occurrence due to piezometric head 
difference between the channel and ground water table. 

Suction/in jectin have also found many applications in modern 
life. For instance, water is sometimes abstracted from lined 
irrigation canals or pipes through perforated boundaries. In 
aerodynamic and internal flows injection has been found effective 


in reducing skin friction and heat transfer to the walls 
(transpiration cooling), whereas suction has been employed to 


reduce flow turbulence and prevent flow separation [Jeroroin, 1970; 
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Fiedler and Fernholz, 1990; Naraaiinha and Sreenivaaan, 1979; 
Subramanya, 1986]. 

In atudyinfi the effect of auction and injection — 
hereinafter auction/ in jection ia used aynonynoualy with uniformly 
distributed suction/injection throu£h the boundary unless 
otherwise specified — two aspects are worthy of note. First, for 
flow over a permeable boundary without any net transfer across 
the boundary referred to as passive porosity — the flow 

properties may be noticeably affected as compared to flow over a 
similar, but impermeable, boundary. This is probably because 

(i) the no-alip condition is not strictly met at the boundary; 

(ii) the no-penetration condition too may be inaccurate for the 
mean velocity (due to possible local variations in U) and 
invalid for the fluctuatin£ velocity component, U’. 

Empirical results indicate that passive porosity increases the i 
skin friction, and chanses the constaints of lo£-law [Beavers and 
Joseph, 1967; Bushnell and hcGinley 1989; Gupta and Paudyal , 1985; 
Schetz and Nerney, 1977; Za£ni and Smith, 1976; Zippe and Graf, 
1983]. Huch of the experimental conclusions reported in ; 

literature are, however, conflictinfi, and some later work eg, 
Richardson and Parr [1991], has indicated that the friction factor 
may increase or decrease with passive porosity depending on 
boundary roughness or pore size. Evidently, a simple 
correspondence between passive porosity and its effect on the flow 
can not be expected. 

i 

Secondly, sel f-preserving flow is not established immediately j 
after commencement of suction/ in jection. Experiments suggest 

1 

that the internal boundary layer that develops from the point of 
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application of injection may not merge with the outer boundary 
layer before a distance ~ 30H to 40H downstream [Simpson, 
1971], Since the establishment of self-similarity with change in 
other wall variables, such as roughness, spans a length of the 
same order of magnitude, ~ 20H to 60H [Smits and Uood , 1985], 
it may be inferred that in the case of suction, as for injection, 
Xd is not less than 20H or 30H. Unfortunately, many experimental 
measurements with suction/injection are at X considerably less 
than 20H eg, Elena [1984] and Fulachier et al [1987]. This limits 
the utility of these data in studying flows with wall mass 
transfer to localized phenomena. 

2.2.2 Fluid dynamic effects of vail auction/injection 

A considerable body of experimental measurements reported 
in literature show that the gross effects of suction/injection in 
TBL flow are similar to those of favourable/adverse pressure 
gradients. These effects are summarized in Table 2.1. Though 
most of the observed effects are similar in TBL and pipe flows, 
some differences have been reported. To cite. Fennel, Eckert and 
Rodi’s [1972] pipe flow experiments with injection showed an 
initial tendency towards laminar izat ion upto X ~ 5H • More 
significantly Eckert and Rodi [refer Narasimha and Sreenivasan, 
1979] indicated that after proper development of pipe flow with 
injection, the akin friction exceeds that of flows without wall 


mass transfer. 



TAble 2.1: Bydrodyn&ai c «ff«c'ta of auction/injection on vail flow 


Parameter 

Effect of suction 

Effect of injection 

^w' 

Increase 

Decrease 

Shear stress profile 

Tends to bulge near 
the wall. 

— 

Mean velocity profile 

Fuller 

More uniformly 
curved. 

Turbulence 

characteristics 

Production, intensity 
and frequency of 
bursts reduce 

— 

Developing length 

— 

~ 30H to 40H 

TBL thickness as 
compared to canonical 
TBL thickness 

Decreases 

Increases 


In an effort to understand the effect of suction on TBL 
flow structure Gad-el-Hak and Blackwelder [1989] applied 
continuous auction and strong, inpulsive suction respectively 
through two wall slots spaced about 100 apart (equal to 

sublayer streak spacing in laminar boundary layer flow) . In the 
first case flow visualization showed a train of hairpin vorti«^ 
behind the slots suggesting a spatial structuring of the near wall 
flow similar to that in TBL flow. Moreover, the boundary layer 
had a fuller velocity profile than the Blasius boundary layer. In 
the second case the main change observed was turbulence burst-like 
events detected by hot film anemometry. In sum, localized slot 
suction indicated some developments similar to turbulent coherent 
structures depending on how the suction was applied. 

There have been several studies on the effect of suction 
and injection on loose boundary open channels, though detailed 
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measuremen'ts of U(y) and T (y) have not been made in most cases. 
Martin [1970] and Martin and Aral [1971] indicated that the 
vertical seepage force, Fg , acting on surface grains of a sand bed 
can be expressed as 

3 P 

“ ^1 Cl + «s> 2.28 

3 Y 

where the coefficient is about 0.35 to 0,5, 6^ is the void 
ratio and 9P / 3 Y is the vertical pressure gradient. 

Inserting Equation 2.28 into Shields parameter Brink and 
Oldenziel [1974] obtained a modified mobility number, , as 



dHp 

a^D [^ - 1 + qci + eg) — -] 

dY 


where a^ is the gravitational acceleration, D is grain diameter 
and Sg = ^ d ^ being the mass densities of grains and 

water respectively, and Hp is piezometric head. The authors did 
not attempt to measure the bed shear stress. Assuming that the 
bed shear stress increases with suction and decreases with 
injection, since the denominator in Equation 2.29 also increases 
with suction and decreases with injection. Equation 2.29 does not 
readily indicate whether sediment mobility is enhanced or reduced 
on applying auction (or injection). However, Brink and 
Oldenziel ’s experimental data consistently show increase in 
sediment transport rate with injection and its decrease with 
suction. An illustrative set of velocity profiles given by the 
authors shows that the velocity profile is fuller for suction and 
more gently curved for injection in comparison to that of uniform 
flow. This agrees, at least quatitatively , with TBL flows. 
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Haclean and Uilletts [1986] conducted OC flow neaaurenents 
with relatively high rates of auction through a porous bed. Their 
indirect measurement of from suction-free grains, indirect 

estimation of from bedload and estimation of from momentum 

balance indicated monotonic increase in T,, with suction rate. 
However, each of these estimates differed widely from the others, 
with the momentum balance estimate being several times the 
estimate from the threshold state of grains shielded from suction. 
The intermediate estimate of was from an equation for bedload 

motion based on sediment mobility similar to Equation 2.29. 
Further velocity measurements by Haclean [1991] for OC and TBL 
flows with high suction rates bring out some basic differences in 
the two flow situations. In contrast to OC flows where |V(y) | 
decreased monotonically with increase in y, in TBL flow V(y) 
bulged slightly in the inner region but was relatively constant 
thereafter . 

Direct measurements of the change in hydrodynamic drag and 
lift forces on surface grains was carried out by Uatters and Rao 
[1971] using 3.75 inch spheres in turbulent flow. Their results 
showed that the drag force decreased considerably with inj-ection 
and increased with suction, while the lift force increased with 
injection and decreased with suction. In both cases the changes 
are more pronounced at low R^, and, at lea^t for the lift force, 
the relative change in magnitude approached a limiting value at 
high R- . Surprisingly, the lift force was found to be always 


negative. 
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2.2.3 Analysia of TBL flow with auction/injection 

Theoretical studies of the effect of wall suction and 
injection on TBL flow £enerally assume the mean vertical velocity 
to be constant (V=Vy) throufthout the boundary layer, and use 
Prandtl’s mixing length method to model the Reynolds stress term 
in the inner region of flow. Host of these models, developed in 
the 1950s and 1960s, have been reviewed in Jeromin [1970] and 
elsewhere. The basic formulation of these methods, assuming no- 
slip velocity at the wall, is outlined below. 

Considering Prandtl’s boundary layer equations to be valid, 

T 3 ” 2 

_ _ = V y.y = + VJJ 2.30 

P 3 y * W 


Though Equation 2.30 can be solved by relating the Reynolds stress 
to mean vorticity through a single equation eg, that of Rotta or 
of van Driest (Equations 2.21 and 2.22), most researchers adopted 
a two-layer model by solving the velocity distributions separately 
for the viscous sublayer and turbulent regions and combining them 
through a common value, Uj^ , at the sublayer-turbulent interface. 
Thus, for the sublayer region Equation 2.30 yields 


tT = -z- [exp(v„ y) - 1] 


and for the turbulent region. 


u 


'*log * 




rs? 

U 


log 


2.31 


2.32 


1 

where uy^^ = (---)ln 
possible changes in 
originally developed 


y + B is the log-law 
the coefficients x 
by Dorrance and Dore, 


velocity admitting 
B. Equation 2.32, 


is the so-called 



2a 


bilofiarithmic law. It shows that for any given y, u is linearly 

CO 

dependent on if x and B are constant. 

Depending on the boundary conditions applied on Equation 
2.32 the coefficients x and B can be determined. Conceptually, 
only one boundary condition is available, which is at the 
sublayer-turbul ent interface (alternatively, at the apparent 
origin of the velocity profile), so that assuming X to be constant 
B can be determined. The various criteria suggested — 
essentially intended for flow over smooth boundaries — are as 
follows [refer Jeromin, 1970; Kocheryzhenkov and Matveyev, 1977] : 

(i) Rubesin [1954] proposed that the Reynolds number, at the 

sublayer-turbulent interface is independent of v^ , 


ei 


UiYi 

y~~ 


= 13.1' 


2.33a 


(ii) Dorrance and Dore [1954] and Stevenson [1963] recommended 

«SJ CN» 


that the origin of u on the wall length scale V /U* ie, y for 

CNJ 

u = 0, is invariant, implying thereby 

B = constant 2.33b 

(iii) Black and Sarnecki [1964] suggested that irrespective of V^, 

fSJ 

uj^ is constant. 


Ui = 11.0 

(iv) Simpson [1966] ignored the modified sublayer 
distribution and prescribed 


2 . 33c 
velocity 


u = y£ = 11.0 * 2.33d 
Dorrance and Dore’s, Rubesin ’s and Simpson’s criteria seem to have 
been more frequently used by subsequent researchers. 
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For the defect re«ion of flow [refer Jeromin, 1970] 
Stevenson proposed u as an extended form of the bilo«arithiiiic law. 
Mickley and Smith suggested a modified wake law by normalizing the 
velocity with the maximum shear velocity in the boundary layer 
rather than the wall shear velocity. HcQuaid used an 
independently computed wall shear velocity taking account of the 
pressure gradient and transpiration velocity. In contrast to 
earlier researchers, Torii, Nishikawi and Hirata did not assume V 
to be constant and, instead related it to U through the continuity 
equation. They subsequently made some assumptions to obtain U(y) 
and T(y). Most of the analytical solutions for the defect region 
cited here show fair to good agreement with empirical data 
[Jeromin, 1970] . 

Kays [1972] studied heat and mass transfer to TBL using a van 
Driest-type mixing length formula for the inner region and a 
mixing length scaling on boundary layer thickness for the outer 
region with a criterion similar to Rubesin’s for the sublayer- 
turbulent interface. The numerical solutions were found to match 
very well with experimental data of velocity and temperature 
profiles including measurements in the developing or relaxing 
length of the inner boundary layer. 

2.2.4 Analysis of OC flow with auction/injection 

Conventionally, the depth average continuity, momentum and 
energy equations are used in the analysis of spatially varied flow 
in open channels. The analysis is valid for inflow to, or outflow 
from, the channel through the top, sides or bottom of the channel 
and at any rate of inflow or outflow so long as the channel flow 



ireina.ins 'turbulen't and 'tha flow doas not chan^a from subcrltical to 
supercritical or vice versa. The essential analysis and main 
results [Chow, 1959; Subramanya, 1986] are discussed below. 

For streamwise increasing discharge, the net momentum flux 
across the channel section is 


+ rASo - ykSf 2.34 

dX dX 

where Fjj and Fp are respectively the momentum and pressure forces, 
y is unit weight of water, A is cross-sectional area, is the 

channel slope (|Sq|<< 1) and S£ is the energy gradient of flow. 
Uith 3P/3X = y k dH/dX and Fj^ = Cj^ Q^/A, where is the 

momentum correction factor and Q is the discharge. Equation 2.34 
gives the water surface gradient as 


So - S£ - 2 C„ Q Q*/(agA^) 
\2 m ^ ^ 3 < 


2.35 


dH 

dX 1 - 0“ T/(agA") 

where Q* = dQ/dX (that is, inflow to channel per unit lengtt^), and 
T is the water surface span. For wide rectangular channels with 
constant Qa the above equation simplifies to 


dH 

So - Sf - 2C„ V 

i '*av Sq 

dX 

1 - 

So 


2.36 

where Vj^ is the inflow rate per unit bed area, and vj^ and u^^ are 
and normalized with ^ agHS^ . Holding S£ and Cj^ constant, 

the water surface elevations can be easily computed if the depth 
and mean velocity at a control section are known. 

For spatially decreasing discharge, it is contended* that 


the energy content of flow per unit mass is not significantly 
affected by water withdrawal. Hence the depth- integrated energy 
equation is applied instead of the momentum equation. 
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dH 



L. 

dX ZagA^ 




dH 

dX 



2.37 


2.38 


for wide rectangular channels with uniform water withdrawal rate, 
where is the energy correction factor. 

The one dimensional analysis outlined above has several 
shortcomings. Firstly, other than H and , many aspects of the 
mean flow and stress field, such as U(y), and T(y), remain 
undetermined. Secondly, there is considerable scope of error in 
assuming Sf based on uniform flow equations, such as Hanning’s 
equation, which may not be valid for gradually or spatially varied 
flow. Thirdly, the coefficients and depend on the shape of 
the velocity profile and may change significantly with 
suction/injection; in fact, is known to be much higher for 
increasing discharges than for uniform flow [Subramanya, 1986] . 
Finally, the analysis does not throw any light on the turbulence 
flow field. 

An attempt to extend the two dimensional analysis of TBL 
flow with suction/injection to OC flows was made by Nakagawa 
[1984], who adopted the bilogarithmic law (Equation 2.32) and 
prescribed a slip velocity at the bed surface. The slip velocity 
boundary condition was estimated by considering velocity 
fluctuations U’ and V’ (produced by pressure fluctuations of the 
turbulent channel flow penetrating the bed) to be correlated. 
The resulting turbulent shear stress in the bed produced 
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streanwise seepage flow, whose magnitude at the bed surface gave 
the slip velocity. It may be pointed out that the author did not 
consider any damping of pressure fluctuations in an uncompacted 
granular medium; nor were possible correlations between U* and U’ 
or between V’ and U' in the bed (which could have induced three 
dimensional seepage flow) taken into account. Most importantly, 
consequent upon Nakagawa’s adoption of the bilogarithmic law, = 
Vy at the free surface. Thus, the free surface must have a slope, 
Sp, relative to the bed. 



This contradicts the author’s assumption of the free surface being 
parallel to the bed. Hence the validity of a self-similar 
solution (Equation 2.32) for the entire flow depth is 
questionable. 

As illustrated by the above discussion there are serious 
difficulties in directly transplanting the transpired boundary 
layer analysis to OC flows subjected to suction/injection. At the 
same time there do not appear to have been any successful two- 
dimensional analysis independent of TBL assumptions. Hence the 
one-dimensional hydraulic equations remain the most effective 
tools for even complex problems involving sucked/injected channel 
flow such as conjunctive modelling of channel flow with subsurface 
seepage or infiltration [Akan and Yen, 1981] and dynamic analysis 
of surface irrigation through furrows [Schmitz and Seus, 1992]. 
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2-3 Hoch&nlcs of lonsi.'tudinAl bodforns 
2.3.1 G«noral 

Turbulent fluid flow with sediment motion over a flat bed 
often leads to a systematic three-dimensional reshaping of the bed 
by locally unmatched erosion and deposition rates until dynamic 
equilibrium is achieved. Among several common sedimentary 
structures formed in such cases, two-dimensional, flow-aligned 
longitudinal bedforms are common to different sedimentary 
bedforms. Some of their observed occurrences in nature are 
[Allen, 19^4; 1982; Karcz , 1967; Hanna, 1969; Kenyon, 1970; 

Glennie, 1987; Werner, et al , 1980] ; 

* in deserts as longitudinal or self dunes; 

* on non-cohesive or weakly cohesive sediment beds in rivers, 
seas and lakes as longitudinal stripes, ridges, lineations, 
(sand) ribbons, or mudgrooves; 

* on bedding surfaces of fine- to medium-grained sandstone as 
parting lineation or primary current lineations. 

Though longitudinal or streamwise bedforms vary widely in 
dimension and, to some extent, morphology, depending on their 
geophysical setting, their common shape and alignment tends to 
suggest that their mechanism of formation and maintenance are 
similar, if not identical. Many researchers have described 
boundary layer as the active agent that sculptures longitudinal 
bedforms on passive, deformable beds, and suggested that 
transverse structuring of the flow, probably in association with 
streamwise helicoidal vortices, are instrumental in this process 
(eg. Allen, 1982; Folk, 1976; Gupta, Parkash and Garde, 1987; 
Hanna, 1969; Muller and Studerus, 1979]. Though many details of 
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the mechaniain that enables a close coupling of the flow and the 
bed may not be understood yet, the importance of longitudinal 
bedforms can be realized from the following facts [Allen 1982; 
Folk, 1976; Glennie, 1987; Tsujimoto, 1989] : 

* Longitudinal bedforms in contrast to transverse (two or 
three- dimensional) bedforms exert little form drag on the 
flow. 

* Knowledge of these structures can be used to infer 
significant aspects of the flow and sedimentological context 
in which they occur. Such knowledge have been used in 
paleocurrent analysis and in identifying desert wind 
direction. 

The above results have accrued mainly from empirical 
knowledge of streamwise bedforms, and further advantages may come 
with better theoretical understanding of their mechanics. As of 
now, two main theories of longitudinal bedforms have been proposed 
by researchers, — one relating them to coherent structures of the 
viscous sublayer, and the other to secondary currents of the bulk 
(or outer) flow. These models are reviewed here. 

2.3.2 Low-speed streak model 

Based on Sorby’s suggestion early this century [refer 
Allen, 1982], parting lineation on parallel-laminated sandstone 
have been conceptually linked with bedload streaks and the 
spanwise organization of quiescent sublayer flow [Allen, 1964, 
1982; Jackson, 1976; Mantz, 1978; Ueedman and Slingerland, 1985; 
Williams and Kemp, 1971]. According to this thesis lateral 
advectlon by the sublayer vortical flow deflects bedload towards 
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low— speed streak regions where, due to their increased 
concentration and low fluid nomentum, the bedload tends to deposit 
in flow-parallel lanes as indicated in Figure 2.1. Jackson [197d] 
and Allen [1982] have suggested that the main impetus for sediment 
motion and the ridge-groove formation comes from turbulent 
bursting events and not the quiescent flow, and that incipient 
bedforms induce lateral fixing of low-speed streak and/or burst 
locations . 


buffer 

region 



The low-speed streak model has been verified in the 
laboratory by Allen [1969], Mantz [1978] and Ueedman and 
Slin£erland [1985]. These studies reveal that the characteristic 
spacing, bedload streaks and sedinent ridges produced on 
the flume bed is about 100 ^ /U* which is the same as that of 
sublayer streak spacings, X Further, shape-fabric studies of 
sandstone deposits shoved that the long axes of grains on parting 
lineations were oriented at small angles to the lineations, which 
is indicative of the impact of weak lateral flow on the lineations 
[Allen, 1982]. 

The evidence cited above does not validate the low-speed 
streak model completely, particularly since the mechanism of 
sublayer flow structuring is insufficiently understood at present. 
Thus, it is yet to be demonstrated that incipient ridge-groove 
formations can control the spanwise organization of sublayer flow 
and turbulence bursts. Without such control, bed lineations will 
disappear almost as soon as they begin to appear without leaving 
any trace. Secondly, it is not very clear why in the experiments 
of Hantz bedload streaks were observed with micaceous flakes 
whereas ripple bedforms developed with granular sediment. 
Thirdly, Gupta et al . [1987] reported lineations with grain shear 

«v> 

Reynolds number, d, in the range 11 < d < 70. This corresponds, 
for a static bed, to transitional bed roughness when the sublayer 
flow structure is at least partially destroyed [Grass, 1971]. 
Hence it is difficult to fully account for Gupta et al’s 
lineations in terms of the sublayer flow. These and other 
<]ue8tions that may arise, such as the possible correspondence 
between the dimensions (length and relief) of lineations and of 


sublayer streaks, underscore the need for detailed hydrodynamic 
and sedimentolofiical studies of lineations and their relation to 
sublayer flow. 

2.3.3. Secondary current model 

Secondary currents in straight open channels usually occur 
as counter-rotating pairs of streamwise vortices whose vertical 
dimensions are on the order of flow depth. They may be caused by a 
number of factors such as asymmetric channel boundary (sidewalls 
or channel banks), the effects of corners and bed disturbances, 
and, in general, any physical condition that produces unbalanced 
hydrodynamic forces in a rectilinear flow configuration. These 
forces result from the inhomogeneity in anisotropic wall 

turbulence and produce Prandtl’s second type of secondary flow, as 
opposed to secondary currents of Prandtl’s first type caused by 
skewing of the mean flow eg., in channel bends [Perkins, 1970]. 
Proponents of the secondary current model for streamwise bedforms 
attribute the existence of the latter to lateral redistribution of 
streamwise momentum and shear stress near the wall by secondary 
flow, with generation of secondary flow being possibly aided by 
the existence of bedforms or lateral grain sorting [Glennie, 1987; 
Hanna, 1969; McLean, 1981; Muller and Studerus, 1979; Nezu, 
Nakagawa and Tomingaga, 1985 ; Tsujimoto, 1989]. 

The secondary current hypothesis is empirically supported 
by the direct or indirect detection of secondary flow cells in the 
presence of lineations, and the correlation of the lateral spacing 
of cells with L, . Flume studies show that the spacing of pairs of 
secondary cells with bed roughness strips sira ulating sand ribbons 
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is equal to * sufieesting that sand ribbons are produced in the 
low nomentuiR regions where bedload tends to deposit [Huller and 
Studerus , 1979; Studerus , 1982 ] . Moreover, sieasurenents by the 
above authors and flow visualization by McLean [1981] showed the 
near— bed flow to be directed away from the roughness strips 
indicating the likelihood of greater bedload concentration and 
lineation development in the smooth zones. 

Nezu et al [1985] reported secondary flow cells in 
experiments with flow-parallel ridges placed on a channel bed. 
They found the near-bed secondary flow to be directed towards the 
ridges which supports the above hypothesis if the ridges are 
considered analogous to lineations. However, the local effect of 
bed ridges should be a decrease in the log-law velocity similar to 
that of roughness strips. Thus Nezu et al ’ s measured direction of 
the near bed secondary velocity is not quite equivalent to those 
of Muller and Studerus and of McLean. 

Unlike the low-speed streak model, attempts have been made 
in the case of the secondary current hypothesis to make 
quantitative analysis or estimates of velocities and shear stress. 
Odgaard [1984] proposed a shear induced secondary flow model that 
makes simplifying assumptions about the distribution of V(y) and 
the profiles of T (y) and U(y). Though these assumptions are not 
rigorously justified and, in fact, give a considerable slip 
velocity for U at the channel bed, the resulting shear stress and 
mean velocity distributions show good agreement with the data of 
Muller and Studerus. 



HcLean [1981] developed a more elaborate matbematical model 
for the mean flow field and stress on a flat bed with spanwise 
periodically varying roughness. Uith an aasumed form of the 
variation of turbulent normal stresses and with an eddy viscosity 
formulation for shear stresses he obtained secondary flow 
velocities as functions of some undetermined coefficients. 
Assuming the secondary velocities and a secondary component of the 
streamwise velocity to be the perturbed mean flow components due 
to variation in bed roughness, he obtained the streamwise velocity 
and shear stress from the X-momentum equation using eddy viscosity 
formula. The unknown coefficients were evaluated by matching 
computed UCy) with his experimental data, and he thus showed the 
mean flow field to be uniquely related to bed roughness variation. 
McLean also suggested that the most likely ribbon spacing is ~ 3H 
conforming to field observations, and that, though the secondary 
circulation is weak compared to the mean flow, the ratio of 
over rough and smooth strips could be as high as 2 or 3. 

In contrast to McLean’s and Odgaard’s approach Tsujimoto 
[1989] assumed secondary flow cells and shoved that it produces 
lateral grain sorting on a non-uniform grain bed. By assuming an 
oblique near-bed velocity vector, and with empirical estimates of 
entrainment rates and saltation distances of different size 
fractions, Tsujimoto’s analysis indicated the development ol 
lateral grain sorting and their attenuation with equilibriun 
transport rates of bedload composed of sand and gravel. 

A 

In concluding, it may be mentioned that whereas the scalet 
of bedload streaks and parting lineations suggest that they ar< 
related to sublayer or wall layer flow structures, seconder: 
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currents o£ the bulk flow seem to be the likely cause o£ sand 
ribbons. However, though secondary currents are purported to 
produce sand ribbons on non-uniform firain beds, such ribbons (or 
lineations or ridges) on graded material beds have not been 
reported in literature. 



CHAPTER 3 


THEORETICAL CONSIDERATIONS AND NODELLING 


3.1 Analysis of flow with suction/injection 

As discussed in the previous chapter the present analytical 
solutions for flows with suction/injection are essentially suited 
for TBL flows. In OC and duct flows the vertical velocity does not 
transpire throuah the boundary layer uniformly ie, V(Y) f V^, and 
the int ermittency of the outer TBL reeion is absent or much less 
in these cases. Herein an attempt is made to analyse spatially 
varied OC flow with uniformly distributed suction/injection and 
with water surface nearly parallel to the channel bed. It is 
assumed that the turbulence generation, diffusion and dissipation 
mechanisms are fundamentally unchan&ed from those of uniform OC 
flow, so that turbulence quantities can be modelled identically in 
both cases. 

3.1.1 Basic formulation 

Consider turbulent flow in a wide, rectangular channel with 
small bed slope, (~ sin 0), as shown in Figure 3.1, subjected 
to injection or suction at a constant rate Vy over sufficient 
length of the channel. For a Cartesian coordinate system with Y 
normal to the bed and X flush with the bed in the streamwise 
direction having origin at the location where dH/dX = 0, let Hq be 
the flow depth and the specific discharge (discharge per unit 
width) at X=0. The following assumptions are made in this 
analysis : 

(i) The flow is incompressible, steady and two-dimensional in the 


mean. 
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Figure 3,1 : Schematic of open channel flow with 
suction or injection 

(li) At X=0 water iminediately below the bed has the same 
longitudinal and lateral pressure gradients as those of the 
overlying fluid, ie, 

3 P 3P 

= = 0 near X=0 3.1 

3 X 3Z 

where P is the pressure difference with gauge pressure. Equation 
3.1 implies constancy of 3P/3 Y across the bed surface, and hence 
constancy of Vy 

Assumming all quantities to be invariant with 2, the 
governing equations of motion for the turbulent flow region in the 
neighbourhood of X=0 are : 

3 U 3V 

3 X 3 Y 


0 


3.2 
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acu^) 

+ 

3X 

3(UV) 

3Y 

3P 9 

+ agSo CU’ 

^ 3 X ® 3X 

9 „ . 





(U’V’ ) 

a Y 

3.3 


3 CUV) 

+ 

3X 

3(7^) 

3y 

3 p 3 — r 

+ ag a;'/) - 

f 3 Y * 3x 

3.4 

9Y 

with the boundary 

conditions 

: 


Ci) 

at Y = Ho , 

V = = 

U’V’ = P = 0 

3.5 

(ii) 

at Y = Yo. 

V = V„ and 

U = 0 

3.6 

where 

ag is acceleration due to gravity and Yq is 

the apparent 


origin of the turbulent velocity distribution. A simple order of 
magnitude analysis indicates that when is small, the equations 

of motion reduce to Prandtl's boundary layer equations albeit, 
with the inclusion of the gravitational acceleration term. 

Herein, some empirical considerations will be used to obtain an 
estimate of Vy for which the boundary layer simplification is 
valid. 

Since dH/dX = 0 at X = 0 and since the bounds on V at Y = 

Yq and Y = Hq are uniform in the neighbourhood of X = 0, it may be 

assumed that 9V/ 3X = 0 for all Y, or V(X,Y) = V(Y1. Hence on 

integrating Equation 3.4 over Y to H^ and differentiating with X, 

1 3P 3V’2 a Ho 3 

- _ / (U'V’) dt 3.7 

P 3X 3X 3Xy3X 

t being a dummy variable. Subsituting the above in Equation 3.3, 



1 1 


a So + ---C-U'V’) + -i-(V'2 

* 3y ax 


+ 


pi **o ^ 

ax Y ax 


= ---(U^D + ---(UV) 

ax ay 


3.8 


From empirical TBL and uniform OC flow data it may be assumed that 
U’ and V" are at most a few times the kinematic shear stress 
which, for uniform OC flow, is agH^S^. Uith water surface slope, 
S = Sq " dH/dX, the third term of Equation 3.8 is thus of order 

d d^H 

-CagHS)] = - a Ho [-^-1 3.9a 

dX X=0 dX^ X=o 



which is numerically small in comparison to the first term if 


Ho d^H 
So dX^ 


<< 1 

X=0 


3.9b 


2 2 

This condition is automatically satisfied if d H/dX = 0 at X=0, 
ie, if the water surface is uniformly parallel to the bed or has a 
point of inflexion at X = 0 . If |d^H/dX^| > 0, it will be shown 
later that the condition is fulfilled if 


I Vy I 

I 1 << 1 3.9c 

I Uav I 

An estimate of the fourth term of Equation 3.8 may be 
obtained by use of Prandtl's mixing length formulation together 
with the assumption of a self-similar X-velocity satisfying the 


continuity equation 
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, - I 3 U I 3 U 

- U’V’ = I I 3.10a 

I 3 Y I 3 Y 

V X 

U(X,Y) = U(0,Y)(1 + ) 3.10b 

Oo 

Subsituting Equations 3.10 in the fourth term and dividing by 

1 3 Hq a 

; C-U’V')dt ~ OC . ] 3.11a 

agSo 3X Y 3X agS^ 


Hence the fourth term may be ignored when 


» V 5 

( « 1 3.11b 


3.1.2. Turbulence stress modelling and governing equations for 
tine-average velocity 

Uith the considerations outlined have. Equations 3.3 and 
3.4 yield 


U’V’) = - a-So + (U^) + — CUV) 3.12 

a Y * 3 X 3 Y 

The Reynolds stress term is modelled here by a gradient transport 

formulation involving kinematic eddy viscosity, 2^^, 

3U 

- U’V’ = 3.13 

3Y 


where V ^ 


|x^y2 (1 


X Y (1 


Y 3 U 

H 3 y 

Y 

— ) U* 

H 


3.14 

3.15 


with the last two equations being the mixing length and parabolic 
eddy viscosity formulae respectively, and Ujt being the local bed 
shear velocity. 

Since 3V/ 3X = 0 and slnc’e HCx) (= H^j) is constant near 
X=0, Equation 3.2 implies that U(X,Y) is a linear function of X. 


Let U(x,y) = F(y) + x G(y) 3.16 

1 

==> V(y) = / G(t) dt 3.17a 

y 

1 

and. / G(t) dt = 3.17b 

Vo 


where x = X/Hq and y ^/^o. Hence from Equations 3.14 or 3.15, 




= li*o + 

Px] 


3.18a 




3F 



where 

C 

o 

= [a; yvT-y - 



3. 18b 




3y 





^ 3G 

3F 

”1 


and 


= L— - / - 

— 

1 ' 

3.18c 



3 y 

3y ■ 




Thus, is a linear function of x. 

On subsituting Equations 3.13, 3.16 and 3.17 in Equation 


3.12, integrating from Y to , and dividing by agH^S^, 


— ~~~~ ^ * * 

^g ^o ^o ^y 


— ) =a'i/) I I + 0fe + X0gg ] 

dy >0 


3.19 


where h = H/Hq , fCy) = F(y)/V^ agH^So , g(y) = G(y) / agH^So , and 
®flf2 following operation on any two function fjCy) and 


that are continuous and integrable in the interval yo^-Y^.^ • 


3.20 


1 1 

Ifl(y) f dt - 2 ; £j^(t)£ 2 (t) dt] 

y y 


Since 2^^ is linear in x. Equation 3.19 implies that hS/S^ should 
be quadratic in x. Thus Equation 3.19 becomes, on assuming the 
following expression for hS/So, where a and a are constants, 


= (1 + ax)^ + a j^x^ 3.21 

^o 

d£ dg 

( + X ) = 

H/agH^So dy dy 

+ ^ O] fl-y) 3,22 

On substituting the turbulent viscosity formula of Equation 3.14 

0 1 2 

or 3.15, and equating terms of order x , X and x respectively, 
three separate equations can be obtained, which, together with 
boundary conditions, are required to yield the functions f(y), 
g(y) and the constants a and *^1 . 

3.1.3 Self-similar solution for U(y) 


(a) Mixing 1 ength formula : Uith Equation 3.14 for t’ Equation 

0 1 2 

3.22 yields the following equations of order x , x and x 



2 2,, 2 

X y f 

= 1 + 0fg 

3 , 23a 


X^y^f ’g’ 

®gfi 

= a+ 

3.23b 



2 


and 

X2y2«.2 

= 

3.23c 


[l t 0£g + x(2at 0gg) 


together with the boundary conditions 
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(i) at y = ,f(y) = «(y) = 0 3.23d 

1 

(ii) / fi(t)dt = Vy 3.23e 

Vo 


where f* = d£/dy, fi’ = dg/dy, v’ = dv/dy, etc. The above 
equations are in integro-diff erential forms, but are easily 
converted to ordinary differential equations, namely, 
d ^ o o 

[X^y^(l-y)f ’^] = 1 + f V + f’CVy- V) 3.24a 

dy 

^22 2 
y Cl-y) f’v”] = - a + [v’‘‘ + v” - v] 3.24b 

dy 

2 


[X^y^Ci-y) v”^] = 3.24c 

dy 

Inspection of Equations 3.23 or Equations 3.24 does not indicate 

2 

any closed form solutions for f(y) and g(y) . However, if v (y) << 

2 

1 then it may be anticipated that 0 £g << 1. Moreover, if gCy) bo 

2 

of order a f(y), then 0 << |a|. Then Equations 3.23 in 

binomial series expansion approximate to 

jCyf’ = 1 + 3.25a 


X ye’ 


d + 


® gg ~ ® f g 


3.25b 


1 

and xve’ = a (1 + 

ft 

Equation 3.25b contradicts Equation 3.25c. The latter is 

2 

as it involves error terms of order x whereas UCy) was 


3.25c 

ignored 

assumed 


linear in x. Making the assumption that the flow is nearly self— 
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aimilar, let 


eCy) = a f(y) 

Thus Equations 3.25a and b imply, 

dg 0. 


3.26 


and 


X y = a + 

dy 

d£ 

Xy = c 

dy 




/ 


3.27a 


3.27b 


Evidently Equation 3.26 is invalid, but, if |0eft| << I al. then 
Equations 3.27a and 3.27b hold approximately. In that case there 
may be error in d(y) as determined from either of these equations, 


but the error will be numerically small as |0 


£& 


« 1 . 


More 


sifini f leant ly, since / g(t)dt is fixed at both extremes y = y^ 

y 

and y = 1, if eCy) be a monotonic function within these limits, 
then the error in determining 0 £g will be relatively small. Thus 
fCy) determined from Equation 3.25a will be fairly insensitive to 
the error in g(y) determined from Equation 3.27a or 3.27b. 

Equation 3.27b with Equations 3.23d and e yield the familiar 
log-law : 


a y 

g(y) = In ( ) 

X Vo 


3.28 


From Equation 3.27b an approximate solution can be obtained by 
taking cue from power laws often used in pipe and OC flows. For 
0 < y < 1 and |nl << 1 the following relation approximately holds, 


/ t‘'dt ~ C )C1 - y) 

y 1 + n 


1 + n'fy 


3.29 



It can 


Ije verified numerically that the error in the above 

equation is ^ ^ 2.2 % for |n| < 0.2. 

Thus it is proposed to seek solutions of the for* 

fl(y) ~ 1 Cy*' - Cq '/■y - cj) 3.30a 

where the constants Cj' and i are determined from the boundary 
conditions 3.236 and d, respectively as 

Cl = yo*^ “ ^o 3.30b 


1 


V 


w 


1+n Vy^ 

(1-yo) c 

1+n 



Cl) 


3 . 30c 


Substituting Equations 3.30 in Equation 3.27a, and equating terms 
of order y^ » y^*^ respectively, the constants n, Cq and a can 

be determined. 


(b) Parabolic eddy viscosity formula : With Equation 3.15 in 
Equation 3-22 and proceeding as earlier, the following equations 
are obtained in place of Equations 3.23a and b : 


df 


n*o 



1 + 0£g 


3.31a 


Thus , 


de 

»y C2c - p) + 0gg - p0fg 

dy 

v/ith relation of Equation 3.26, 


3.31b 


gCy) =CCf(y), 


u*o 


3dy "■ 


dg 

dy 


a + 0gg 


3 . 32a 
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and tin^Xy a 

dy 


■^L L. -i'vMKlf 
, 1. T., KA^;PUR 

im.lio.Pi.imU 


3.32b 


whose solutions can be obtained in the form of Equation 3.30a as 
in the previous case. 

(c) Velocity Profile ; The solved dimensionless velocity 
distributions for the mixing length and parabolic eddy viscosity 
models can be represented in the common form : 


«(y) 


u*o 


In (— ) 


fiCy) = 1 [y*^-yo^“6n^ C '^yn)] 


3.33 


3.34 


fCy) = in ty*'-y^-6n^C '^y- '^yo)3 


3.35 


where 


Vy [l-n-Sn^-c^] 

a^XviHo Cl-yo ) [ l-n-3n^ +n Vy^, -ci ] 


Co ~ 

Cl = yo*' - Co '^yo 


(l-Vo) Cl-n-3n^+n 


3.36 


with a = 


aCl-yo) (l-n-3n^+n ^Vq-Ci) 
2(l-n-3n*-ci) (l-n-3n*-0 .5ci ) 


2 for mixing length formula 
1 for parabolic eddy viscosity 


It can be verified numerically that Equation 3.34 approaches 
Equation 3.33 and Equation 3.35 approaches the law of the wall in 


the limit as n — > 0 or, as v„ — > 0. The likely error in 


estimating g(y) can be reduced by taking the average of Equations 
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3.33 and 3.34. Alternatively, ft(y) can be computed by replacinc 
the exponent n by nj^ = n/2 in Equation 3.34 and evaluating the co- 
efficients cq, cj and 1 accordingly. 


(d) Parameters a and : Since there is an error in the assumed 
equality gCv) = ttf(y). the parameters a and (3 representing water 
surface curvature and rate of change of bed shear velocity, the 
parameters a and P will also absorb the error. To reduce the 
scope of error in a , the depth-average value of g(y) and f (y) 
may be used in Equation 3.26, so that 

1 

/y «Ct)dt v„ 

a = = 3 3, 

/ f(t)dt U^y 

dg/dy 

Substituting P = [ ] y_yQ in Equations 3.25 or in Equations 

df/dy 

3.31, P is found to be approximately, 


2.(a-l)0^ 

3/a~ [ «l-_] 

2^^fg,o 

where is the value of 0fg at y = y^. 


0 


fe,o 


~ -2 v„ Uav 


3 . 38a 


3 -38b 


(e) Comments on the analysis ; CiD Uhen using the mixing length 

2 2 

model , terms of order 0 £g, 0 gg and higher were ignored in 

Equations 3.25. This may give rise to significant error in the 

2 

solution set unless 0 gg << 1. The solution set for the parabolic 

A 

eddy viscosity formulation is free of this error. 

(ii) From Equations 3.21 and 3.37 the water surface curvature at 
x=0 is found to be positive for suction and negative for 


injection. 
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d\ 1 

-2 Vy 

So 



___ , 


— 


3.39 

dx'^ 1 x= 

0 «av 




From the above 

relation it 

follows that 

the 

inequalities of 

Equations 3 . 9b 

and 3.11b 

necessary for 

the 

boundary layer 


approximation are fulfilled when IVy/u^^^l << 1. 

3.1.4 Evaluation of 

The velocity distribution derived above is incomplete 
without defining the apparent origin of the turbulent velocity 
distribution. Towards this end a two layer model is invoked for 
hydraulically smooth bed flow assuming uniform transpiration rate 
in the sublayer and an identical mean velocity, u^, for the 
sublayer and turbulent region flows at their interface, y^ . Thus, 
for the sublayer region 3u/ 3 x =0, and momentum balance in the 
x-direction gives 


2 

R* (u + VyU - y) 

dy 


3.4 0a 


where R. 


V 


3.4 0b 


is a Reynolds number. On integrating Equation 3.40a over y, and 
letting u = Uj^ at y = yj^ , 


1 [exp CR*Vyy) - 1] y 
u = ) + 




R*Vy V^ 

[Vy U^ - y^] 

[exp (R* Vy y^) - 1] 


3.41 


3.42 


R* V 


w 
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An independant relation between Uj and may be defined by one of 
several sugeestions made by researchers cited in Chapter 2. 
Herein Rubesin’s criterion is adopted with a modified value at the 
turbulent-laminar interface corresponding to the values of X and B 
in the log-law in a given uniform flow situation. For example, 
with X =0.4 and B = 5.5, 

U£ y^ — 135.373 

135.373 

or , ^1^1” 3.43 

R* 

Thus, if Ut (2 is known then Equations 3.4 2 and 3.43 together yield 
the values y^ and Uj_ from which y^ can be determined ~using 
Equations 3.35 and 3.36. Since u*q is also implicity defined by 
Equation 3.35 and 3.36, can be determined from the above 

equations. However, since u*^, y^, n and Cj^, are interdependent, 
an iterative scheme has to be used for determining the various 
constants. A suitable procedure is outlined below. 


Step 1 


Step 2 


Step 3 


Step 4 


Step 5 


: Given X , B, R* and v^, compute initial values of y^ 
and y£ from uniform flow criteria. 

: Compute u*^ from Equation 3.42 and U£ from Equation 
3.43. 

: Compute n, c^ and c from Equations 3.36 and 3.34 with 
f Cyx ) = . 

: Compute y^ and from Equations 3.36 and 3.48 (given 

later). 

: If y^ and Muq almost same as in Step 2 then 

terminate computation; else repeat Steps 2 through 5 
with new value of y^ . 
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It may be noted here that the shear stress distribution 
derived from Equation 3.40a for the viscous sublayer may be 
different from that of the turbulent region, so that the shear 
stress of the two regions at the interface is not equal. This 
apparently contradicts the notion of a two layer flow with common 
properties at the interface. However, the two layer approach is 
essentially an approximation to overcome the problem of exact 
stress modelling near the bed. The discontinuity in stress can be 
overcome by using a van Driest type of mixing length formula 
for example, Kays, 1972], but this is not attempted here in order 
to eliminate analytical difficulties. As will be seen in Chapter 
5 the divergence between "laminar” and turbulent shear stresses at 
the interface, assuming identical bed shear stress for the two 
cases, is not sufficiently large to warrant a modified eddy 
viscosity formulation. 

For hydraulically rough flow regime (kg > 70) the sublayer 
flow is believed to be completely destroyed or, at least 
inconsequential in defining the velocity distribution origin. 
Hence it may be assumed that, for a given bed roughness y^ is 
independent of v^ . For the transition regime a priori definition 
of y^ in the absence of empirical evidence is difficult. 

3.1.5 Computation of flow parameters 

(a) Ilean velocity and flow depth ; From Equations 3.35 and 3.36 

the depth-average velocity, , is found to be : 

2 

Uav ~®*Cl~n-3n-c^) 3.44 

Thus, given the specific discharge at x=0, the flow depth is 
easily obtained. 
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Cb) Vertical velocity : The distribution of the normal velocity 
component over y is obtained by integrating alternatively Equation 
3.33 or 3.34, 

fv^^[(y-l) (1 + Iny^) - y Iny] 


vCy) 


[y^ - 1 - Iny^,] 

Vy (l-y) (1 - n - 3n^ + n '^y-cj^ ) 
Cl-yo) (1 - n - 3n^ + n 


3.45 


3.46 


(c) Shear stress distribution : From the velocity distribution 
and the eddy viscosity model used, the shear stress distribution 
is derived. 


T'*‘(y) = (l-y) [ran3C(y*'*-3n N/’y) 1^ 3.47 

===> u*o ~ mnKCyQ*' - 3n Vy^) 3.48 

where = T/( . 


(d) Uater Surface Profile ; From Equation 3.21, 

dh 

(Sq ) - So (1+ ax)** 

dx 


3.49 


Using P icard integration and ignoring terms of order x and more, 

M X 


h ~ 1 - a [1 + 


■] 


3.50 


3 3 

which may be considered accurate for (ax) and (SqX) numerically 
much less than 1. 


3.1.6 Extension of results 

It may be noted that while the present analysis gives an 
estimate of the hydrodynamic drag force on granular bed material 
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in a loose-boundary channel, it does not help in estisiatin£ the 
hydrodynamic lift force acting on bed grains. The lift force, 
which plays an important role in bed grain motion, may not be 
significantly affected by changes in the channel flow. Rather, it 
should depend on the vertical seepage pressure gradient as 
suggested by some researchers [eg. Brink and Oldenziel, 1974 ^ 
MacleaiT and Willetts, 1986]. 

The present analytical results should hold good also for 
duct flows, though with reservations about the applicability of 
simple eddy viscosity models for the core flow region. In duct 
flows the question of free surface curvature does not arise, and 
V = 0 at the centreline for all X. Thus, for flow in rectangular 
ducts with uniform suction/injection through opposite boundaries, 
and with 9 P/ (Yd X) replacing Sq, the above results should hold 
identically. Since the flow is only approximately self-similar, 
the parameters Vy , u*^ , CX , 3, and the functions f(y) and g(y) 
would need to be evaluated continually (or at intervals) over X. 

In many field situations involving seepage or infiltration 
through channel bed, the seepage rate may not be precisely 
constant. For the case of slowly varying suction or injection in 
open channel flow, dVy/dx = X , 

1 

v(x,y) = Vy (1 + X x) fgCt) dt 3.51 

y 

Xx 

and u(x,y) = fCy) +xe(y) + ^ 

2 

The governing equations may be rederived as earlier to give 


u*Q/C.yg’ = 2a- 3 + 0 (g+ xf)e ■ ^ 


, a=l 


3.53 
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Further approximations made herein are : 

(iv) The flow is of high aspect ratio, and is incompressible, 
uniform and steady in the mean. 

(V) Turbulent stresses and can be modelled in terms of 

an eddy viscosity and mean vorticity components in the z and 
y directions, respectively. 

("vi) The mean streamwise and vertical velocity components vary 
periodically with Z in consonance with the variation of bed 
roughness . 

The analytical approach here is similar to that of HcLean 
[1981], but no attempt is made here for an exact solution of the 
cross-stream velocity components. 

3.2.2 Formulation 

Consider flow over a ribbon bed with mean depth, H^, as 
indicated in Figure 3.2. The shape, strength and sense of the 
secondary flow cells are unknown. For the simple case of a 
sinusoidal variation of In y^ , where y^ is the dimensionless 
apparent origin of velocity distribution above bed level. 


In y^ = In y^^ [1 + e cos (kz)] 3.55 

where is the Z-averaged value of In y^, k = 27T /Zo is wave 

number with Zq as the wavelength of ribbon spacing, and z = Z^/H. 
It follows from assumption vi of the last section, 

U(y,z) = U^Cy) + U^Cy) cos (kz) 3.56 

and V(y,z) = V^Cy) cos (kz) 3.57 

Considering all quantities invariant in the X-direction. th 
equations of motion are : 




61 


where velocities have been normalized with , lengths with 

Hq , and pressure (P), normal stresses (O) and shear stress ("C ) 
with rH^So, 

p = P/C/HoSo), (T* = a/(/HoSo),T+ =T/(yH^S^) 

Then from Equations 3.57 and 3.58, 


- 1 dVo 

w = sin(kz) 3.62 

k dy 

For the above set of equations, invoking a parabolic eddy 
viscosity formulation for modelling the turbulent stresses and 

namely. 


3u 

"’•xy ■ u*;6y(l-y) 

3y 


3.63 


3u 

■^^xz ~ u* xy (1-y) 

3z 


U* 

implies u* 1 + u*j^ cos(kz) 


X y om (^“yom^ 



y om 


= 1 


3.64 


3.65a 


3.65b 


duj^ 

X y om ^^~yom^ ~ 

dy 


y om 


= u*i 


3.65c 


Hence an order of magnitude analysis of Equations 3.58 and 3.59 


indicates 
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til (y) 
tio(y) 


3.66 


V ^ w ^ E 


If y , and *^yz ar« also modelled similarly as *^xy and "^xz 
in terms of an eddy viscosity and the mean velocity gradients 
3 v/ 3y, 3w/ 3z and 3v/ 3z + 3w/ dy, then Equations 3,60 and 
3.61 indicate that 


3.67 


Hence the lateral variation of flow depth, AH, normalized by Hq, 


is of order 


which is insignificant for small bed slopes. 


3.68 


3.2.3 e - order solution 

On substituting Equations 3.56, 3.57 and 3.62 through 3.65 
in Equation 3.59 and simplifying. 


d au^j a au^ 

— [jcy (1-y) 3 + iXY (1-y) u*i ] 

dy dy dy dy 


dUo dui 


+ [ o; y (1-y) (ti*i + )3 cosCkz) 

dy dy dy 


d aui 2 

+ [JCy Cl-y) u*i ]cos (kz) = 

dy dy 



63 



- 1 + 

1^0 

dUo 

k^a:y Cl 

-y) 

>^1] 

cosCkz) 







dy 







4- 

2 Ik^X 

1 y 

Cl-y) u 

*1 - u 

1 

-!2_ 

] cos^(kz) 


3. 

69 







dy 





Separating terms o 

if order 

cos(kz) , 

1 2 
cos Ckz) and cos 

(kz). 



d 



du 

o 

d 



du 




— [ 

X y (1- 

y) 

— ] 

+ [ X 

y(l 

-y) 

u»i ] = 

1 

3, 

. 70 

dy 



dy 

dy 



dy 




order : 

[1] 




[ ! 


[ 

1] 



■,d 



dui 

d'lo 







— 

C yCl- 

■y)( + 

u*i 

•)] 


, --- f 

yCl- 

•yDui 


dy 



dy 

dy 



dy 







t e] 




[ 6 ] 

te 

] 3 

.71 


[XyCl-y)u*i ] = 2k^ p( yCi-y)u*i - Zu^ 3.72 

dy dy dy 

order : te^] [ 6^ ] [ e* ] 

On comparing the orders of magnitude of each term of the above 

2 

equations as indicated, if e << 1 then Equation 3.72 can be 
ignored altogether, and the second term of Equation 3.70 also 
vanishes. Hence the latter yields the conventional log-law. 


1 y 

yom 

and Eq .uation 3.71 simplifies to 


d du^ Vq 

— IXy (1-y) 1 - k xy d-y) u*i + — 

dy dy ^ y 


3 . 73 


3 . 74 


The boundary conditions for the above are 
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Ci) 

at y = y^, u^Cy^) 

e 

3.75 



du. 


(ii) 

at y = 1 , y(l-y) - 

-i = 0 

3.76 



dy 



From the above set of equations, if v^fy) is known then 
ujCy) can be uniquely determined. The determination of v^fy) is 
however dependent on the modelline of turbulent stress terms of 
Equations 3.60 and 3.61. The conceptual difficulties attendant 
upon such modelline necessitate additional assumptions such as 
eddy viscosity formulae for these terms [cf. McLean, 1981]. Since 
the order of maenitude of these terms is small, errors in the 
modelline of these stresses in terms of the mean flow field are, 
however, likely to be significant. Partly for this reason and 
partly as computational shortcut, it is proposed to assume the 
form of VqCy) as 

Vo(y) = ao y®(l-y)“ 3.77 

where m > 1, n > 0 satisfying the boundary conditions 

dVo 

(i) at y = y^ , v^fy) = = 0 3.78 

dy 

(ii) at y = 1 , v^Cy) = 0 3.79 

Suitable values of the exponents m and n may be assumed based on 
empirical data. To determine the coefficient a^ of Equation 3.77, 
a second boundary condition on uj^(y) may be defined based on 
empirical data, eg. , 

ui(y) = C uiCy^) 

where C is a constant of proportionality. 



3.2.4 Limltfiitidna 

The above analyaia ia an approximate one. It ia eaaily 
seen that, if e la not very small, then 

(i) Bquation 3.70 shows that u^Cy) does not follow the lofi"law; 
(11) liquations 3.71 and 3.72 may not ftive Identical aolutlona for 
hl(y)- 

It follows from (11) that the simple form of uCy,a), which**'*'‘hai8 
only one sinusoidally varylnfi component of z, ia not valid. From 
a practical point of view lateral bed level variations due to 
projection of ribbons above mean bed level will also affect ^ the 
secondary flow. Iforeover, the assumption that In y^ varies 
sinusoidally with z may be quite erroneous in the tt^suc-al case. 
Taking account of these factors, refinements in the above analysis 
may bs neither useful nor easy In the absence of adequate 
experimental data to verify all analytical premises. 

In Chapter 5 the analytical model has been calibrated with 
experimental data of Muller and Studerus [l979] in order to estim- 
ate the bed shear stress variation due to sand ribbons. The re^- 
distributed siear stress as computed has been employed hence in 
hydraulic and sedimoritologicai analysis of sand ribbons generated 
in experiments described in the following chapter. 



CHAPTER 4 


EXPERIHENTS 


To test the hypotheses and seek ansvers to the queries 
outlined in the previous chapter, two sets o£ experiments were 
conducted, one set with immobile bed to test the analytical 

model for suction/in jection, and the other set to examine the 
morphological characteristics and hydrodynamic implications of 
sand ribbons and ridges. The experiments were carried out in a 
glass-sided tilting flume 8.0 m long, 0.4 m wide and 0.6 m high 
and provided with a 0.2 m long approach section. A honeycomb 
structure was used at the inlet to provide rectilinear flow, and a 
0.7 m X 0.4 m flat plate with 3-4 mm gravel glued on it was placed 
at bed level at the inlet to ensure developed flow in the main 
flume section. The discharge in the flume, supplied from an 
overhead tank, was regulated by inlet valves and measured by a 
calibrated V-notch in a tank downstream of the flume. The flume 
water level was adjusted by tailgate and tail wier plates, and the 
water and bed levels were measured by point gauge of least count 
0.1 mm. All experiments were conducted at subcritical flow 
conditions . 

4.1 Experiments with suction/injection 

The main aim of these experiments was to find out the flow 
depths with uniformly distributed suction/injection in a region 
where the water surface profile is approximately parallel to the 
bed. For this purpose a fine wire mesh was spread on “a coarser 
screen covering a 7.0 m long and 0.15 m high structural frame 
placed on the floor of the flume. The two ends of the frame were 
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sealed with flat plates spanning the flume width. An inverted 
filter, 15mm thick, with coarse sand to medium gravel was laid on 
the wire mesh. A sand bed with medium sand size D 50 = 0.53 mm 
and average thickness of 150 ram was prepared on the filter. 
During experiments on OC flow with suction/in j ection , the sand 
layer was topped with two to three layers of fine gravel, 2—4 mm 
diameter , which allowed for higher critical shear and normal 
stresses on the bed than the sand surface. Figure 4.1 shows 
relevant details of the experimental set-up. 

To apply suction/in jection to the channel flow, water was 
abstracted/supplied at a constant rate through three 76 mm 
diameter inlets at the flume bottom. The withdrawal /supply rate 
was effected by an open tank. Tank A, and adjusted by a common 
valve for all three inlets (see Figure 4.1). In experiments with 
injection water was pumped to the tank at a constant rate, 
adjusted to maintain the desired supply to the flune. For runs 
with suction, water was sucked from the flume by operating an 
outlet valve from the tank and measuring the volunetric outflow 
rate. In both cases the net inf luent/ef fluent seepage rate was 
cross-checked with the difference in measured flume outflow rates 
with and without seepage. 

The water supplied by Tank A produced a consxant pressure 
reservoir in the bottom 150 mm height of the flume ; excluding the 
volume occupied by the frame). Since there will be a streamwise 
head drop in the channel flow due to free surface gradient 
(dH/dX), the seepage pressure head across the sand bed may be 
expected to vary with X. As mentioned by Pennell et al [1972] for 
pipe flows, this may result in non-uniform sucti on/ in j ection rate 
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in the X-direction. To ensure uniform the sand bed thickness, 

, was varied in order to provide uniform vertical seepage 
gradient, 9 Y, in every run. Towards this end, first 
suction/injection was applied to a uniformly thick 150 mm), 
saturated sand bed with almost standing water in the flume. 
Pressure tappings in the flume walls at 1100 mm X-intervals, 30 mm 
Y-intervals and penetrating about 30 mm into the sand layer, as 
shown in Figure 4,1, indicated that 9Hp/ 9y was independent of X, 
The mean V^, computed from the net inf luent /ef fluent seepage flow 
and the eand surface area (6300 x 400 ) , was found to follow 
Darcy’s law. 



^Hp 


il 

> 

- K — -- 

4.1a 

3Y 


K = 

0.268 0.011 mm/s 

4.1b 


for - 1.05< 9Hp/ 3Y<1.20 where Hp is the local pressure head, and 
K is hydraulic conductivity of the sand bed. For 9Hp/ 3Y > 1.20 

was less than that predicted by Equation 4.1 a sufigesting 
nonlinear effects (non-Darcy flow). For 3 Hp /3 Y <-1.05 the bed 
surface become unstable with boiling of sand at sporadic 
locations. Hence, to maximize suction/injection rates within 
practical limits, the seepage pressure gradients used were 
9 Hp/9 Y ~ - 0.80 corresponding to V„ ~ 0.215 mm/s for injection, 
and 3 Hp/ 3 Y ~ 1.30 corresponding to = -0.345 ram/s for suction. 

Uith application of suction/injection there was slight 
contraction/expansion of the sand layer. The sand bed thickness, 
Htj(mo^waa found to vary with the seepage pressure gradient 
approximately as. 




Figure 4,1 i Schematic of experimental set»up 
(a) cross section; 
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Hv 


9K^ 9H_ 

- 1.8 (1 - 0.2 -) 

9Y 3Y 


4.2 


for -1 ^3Hp/3Y^1.6 where is the initial bed thickness. 
The variation of bed thickness at constant 3Hp/3 Y was fairly 
uniform over the bed. 

For each experimental run, assumine dH/dX of OC flow equal 
■to ~ Sq , the sand bed thickness, Hjj(x), required to maintain 
constant 3Hp/9 Y along X was computed for = 150 mm at X = 0 , 


dY 


Hb = Ht,o - S^X 

dHp 


or, Htj 


r 150 + 1.25 SqX 


: injection 


4.3 


4 . 4a 


^150 - 0.77 SqX ; suction 4.4b 

Before every run the sand bed thickness was adjusted to that of 
Equations 4.4. With quiescent flow in the flume, several cycles 
of injection and suction were applied, each cycle lasting about 15 
to 60 minutes, in order to expel entrapped air from the sand bed. 
The process also produced slight compaction of the grain matrix. 
The sand permeability and bed expansion/contraction relations of 
Equations 4.1 and 4.2 are for this type of equilibrium bed 
condition. 

On acheiving stable conditions, a few layers of 2-4 mm 
gravel were laid on the sand bed. The bed surface was smoothened 
with a wooden flat and bed levels measured at section LI to L8 
indicated in Figure 4.1. Flume discharge from the overhead tank 
was adjusted to the desired initial value, and suction/ in j ection 
was applied. The tail gate was raised or lowered to vary the 
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flume water level. Uith every chanee in tail fiate elevation the 
water surface profile, the seepage pressure gradient and the water 
level in Tank A took about 10 to 30 minutes to adjust to a 
mutually stable configuration. It was observed that while the 
water surface could have a rising, falling, or mixed profile, the 
profile was usually curved. When parallel to the bed near the 
middle of the flume the water surface had convex curvature (convex 
upwards) for injection and concave curvature for suction. In 

order, to determine the self-preserving flow parameters in the 
region where dH/dX = - S^, the location of optimum water level 
(the origin X=0 in Figure 4.1b) was chosen at 3575 mm from the 
point of commencement of suet ion/ in j ect ion. Assuming the 

developing length of the flow with suction/ in jection to be less 
than 40H (vide Chapter 2), since H < 35 mm in all runs, the flow 
in the region - 2000 < X < 2000 mm could be assumed to be nearly 
self-similar. 

In most runs, the water surface curvature was not 

sufficiently high as to determine the location of optimal flow 
depth precisely. Hence, a set of profiles with dH/dX X=0 

were recorded by slightly adjusting the tailgate position, and the 
profile with dH/dX = 0 at X=0 was interpolated from these data. 
The outflow rate from the flume was also measured and its 
difference with the discharge at flume inlet, Q^, was checked with 
the inflow/outflow rate from/to Tank A. From the total sucked or 
injected flow rate, Q, the suction/injection rate, V„, 
flume discharge at X=0 were computed. 


and the 
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6300 X 400 
3575 


Qo = Qi + Q 

6300 

where Q is the discharge per unit width of flume in ma^/s. Bod 
levels with suction/injection were measured at all sections H 
through L8 and verified against Equation 4,2. Finally, the 
discharge Q^y was run through the flume, and the normal flow depth, 
was determined. The water temperature was recorded during 
every- experimental run. 

All relevant experimental data are recorded in Tables 4.1 
and 4.2. The table also gives the equivalent grain roughness. Kg, 
for the gravel bed computed for uniform flow Qq from the log-low 
assuming X - 0.40 and B = 8.5 (hydraulically rough bed). 

4.2 Experiments on longitudinal bedforms 

The experiments described here were conducted in order to 
determine the conditions of existence, morphological 
characteristics and hydrodynamic implications of ribbons or 
lineations on sand bed with equilibrium sediment transport. The 
experimental set-up is the same as that of the previous section, 
save that the sand bed composition was different and that, the bed 
was not overlaid with gravel or other material. Three different 
compositions of natural river sand, designated SI, S2 and S3, were 
used in these experiments. Whereas SI was essentially fine to 
medium sand, S2 and S3 had increasingly higher amounts of coarse 
sand, with a small amount of very fine gravel being present in S3 
sand. Thus both the medidlk diameter and grading of the sands were 
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Table 4.1 : Observed Open Channel Flow Variables With 
Influent /Effluent Seepage 


Run 

Si 

mm /s 

Vv 

mm/s 

So 

mm /s 

So 

XIOOO 

»n 

tm 

Ho 

inin 

V 

cp 

Ks 

mm 

R1 

7097 

-0.339 

5885 

4.26 

21.4 

20.65+0.2 

0.93 

6.02 

R2 

7928 

-0.346 

6691 

3.80 

23.3 

22.05+0.2 

0.93 

5.25 

R3 

7160 

-0.330 

5980 

2.82 

24.1 

22.9 +0.3 

0.92 

5.73 

R4 

8782 

-0.335 

7584 

2.36 

28.7 

26.95+0.3 

0.92 

5.27 

R5 

6491 

-0.352 

5233 

1.76 

24.9 

23.1 +0.5 

0.91 

4.79 

R6 

6342 

-0.340 

5127 

1.22 

27.2 

25.0 +0.5 

0.91 

4.62 

R7 

5845 

0.208 

6589 

4.40 

22.5 

23.2 +0.3 

0.93 

5.84 

R8 

6389 

0.224 

7190 

3.72 

24.5 

25.15+0.3 

0.93 

5.37 

R9 

' 5453 

0.216 

6225 

2.80 

24.3 

25.1 +0.3 

0.93 

5.37 

RIO 

5470 

0.208 

6214 

2.34 

25.4 

26.5 +0.4 

0.91 

4.90 

Rll 

4653 

0.219 

5436 

1.68 

25.6 

26.6 +0.4 

o:9o 

4.54 

R12 

4740 

0.210 

5491 

1.20 

28.7 

30.1 +0.6 

0.90 

4.95 


Table 

4.2 : 

Observed Water Surface 
Inf luent /Effluent Seepa£e 

Elevation (in mm) 

: Optimum Elevation at 

With 

X=0. 

Stn . 

XCmm) 

Run 

LI 

-3575 

L2 

-3125 

L3 

-2500 

L4 

-1875 

L5 

-625 

L6 

0 

L7 

625 

L8 

1875 

L9 

2500 

R1 

23.6 

23.0 

22.2 

21.6 

21 .0 

20.65 

20.8 

21.4 

21.9 

R2 

25.2 

24.2 

23.2 

22.6 

22.1 

22.05 

22.1 

22.6 

23.7 

R3 

25.9 

25.0 

24.0 

23.4 

23.0 

22.9 

23.0 

23.4 

23.9 

R4 

30.4 

28.8 

27.9 

27.4 

27.0 

26.95 

27.0 

27.3 

27.7 

R5 

26.6 

25.2 

24.1 

23.5 

23.1 

23.1 

23.1 

23.5 

23.8 

R6 

30.1 

27.6 

25.8 

25.3 

25.0 

25.0 

25.0 

25.3 

25.6 

R7 

21.3 

21.8 

22.3 

22.8 

23.2 

23.2 

23.2 

22.8 

22.4 

R8 

23.3 

23.9 

24.5 

24.8 

25.1 

25.15 

25.1 

24.9 

24.5 

R9 

23.1 

23.8 

24.4 

24.7 

25.0 

25.1 

25.05 

24.7 

24.5 

RIO 

24.5 

25.2 

26.0 

26.2 

26.5 

26.5 

26.5 

26.25 

26.0 

Rll 

24.0 

24.6 

25.1 

25.3 

25.6 

25.6 

25.6 

25.3 

25.1 

R12 

26.4 

27.4 

28.1 

28.5 

28.7 

28.7 

28.7 

28.5 

28. 4 
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increased from SI to S2 to S3. The £rain size distribution of 
SI, S2 and S3 sand samples based on sieve analysis are shown in 
Figure 4.2. For each sand type Kramer's uniformity coefficient, 
Cj^, shown in Table 4.3 was computed from the formula [Garde and 
Raju, 1985] : 

50 

Ck = 4.5 

100 

E D£ Afi 
50 


where Af ^ is the percentage interval on the ordinate of Figure 
4.2 having representative (mean) size, Dj. 


Table 4.3 Summary Data on Sand Bed Composition 


Desig- 

nation 

Sand Type 

Size 

range 

(mm) 

Di6 

(mm) 

°50 

(mm) 

^84 

(mm) 

*^90 

(mm) 

Ck 

SI 

F ine 

to medium 

0.09 to 

0.15 

0.28 

0.445 

0.50 

0.48 


sand 


0.70 






S2 

Fine 

to coarse 

0.09 to 

0.255 

0.52 

1.05 

1.31 

0.35 


sand 


1.00 






S3 

Fine 

sand to 

0.09 to 

0.23 

1.14 

2.38 

2.50 

0.18 


very 

fine 

3.00 







gravel 


In each of the present experiments with suction/injection, 
the variation in sand bed thickness was partly carried out in the 
early part of the run by adding or removing sand at the downstream 
side while raising or lowering a plate embedded spanwise at the 
end of the sand bed. With flume discharge at about incipient 


finer 


5 
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condition, some amount of suspended material was initially washed 
away from the bed surface. Experiments with SI sands were 
conducted at just above threshold conditions. It was observed 
that any significant sediment motion produced ripple bedforms, and 
neither suction nor injection promoted the development of any well 
developed lineations. In many experiments, however, fine streaky 
structures of very low relief were perceptible on the sand bed. 
These were in the form of longitudinal grooves, a few centimeters 
in length and of very small transverse and vertical dimensions (on 
the order of grain diameter). These striations often occurred in 
clusters on the bed, sometimes bifurcating or merging with 
neighbouring striations. Their visibility was slightly improved 
by adding a little chalk dust to the flow just before terminating 
the run so that the chalk powder tended to settle into the 
grooves. The only characteristic dimension of these striations 
that could be meaningfully measured was their transverse spacing, 

The experiments with S2 and S3 sands were conducted with 
stable sediment transport rates. The bedload at the downstream 
sink was collected at every 15 minute intervals and fed gradually 
at the upstream end. During the first hour of the run the water 
and bed levels had to be adjusted slightly in order to maintain 
(nearly) uniform flow depth and seepage pressure gradient when 
suction/ in j ection was applied. 

In experiments with S2 and S3 sands, lineations or ribbons 
developed on the sand bed at modest bedload transport rates. The 
lineations, which commenced at or near the upstream end, ran down 
the entire flume length, being usually more well developed in the 
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downstream two-third length of the flume bed. The lineations 
appeared consistently uniform in size, spacing and surface texture 
after a few hours of run. They were characterized by both relief 
(ridge - groove combinations) and grain sorting (coarse and fine 
stripes). Figure 4.3 shows a typical photograph of sand ribbons 
developed on the flume bed. Though the applied suction/injection 
rates were less than those used in Section 4.1, did have some 
perceptible influence on both sediment motion and bedform 
characteristics. In general, injection increased bedload while 
suction had the opposite effects. For S2 sand the longitudinal 
bedforms tended to disintegrate into primary ripples as described 
by Hantz (1978) at high bedload transport rates. Visually it was 
observed that moving bed grains mostly adhered to the sloping 
ridges rather than the valleys between ridges. 

All experiments were terminated within 4 to 6 hours after 
stable flow, bedload and bedform conditions appeared to have been 
achieved. The flume water surface profile, water temperature and 
bed measurements were carried out in all runs. The latter 
included the lateral spacing between lineations (L^) and the 
elevations of crests and troughs of the ridge-groove formations. 
The mean difference between the elevation of a crest and its two 
neighbouring valleys was designated as crest height, Ly . Sediment 
samples scraped from the surface of ridges and grooves 
respectively, and those sampled from the bedload collecting in the 
sink were oven-dried, and their sieve analysis was carried out. 
The experimental results are given in Tables 4.4, 4.5 and 4.6. 
For SI sands the bedload was negligible. In other experiments the 
bedload (dry weight) transport rate, Qg, is given in kg/ro/hour. 
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Tabl« 4.4 : Linaations on SI sand bed 


Run 

(inni^/s ) 

Ho 

Cmm) 

Cep) 

So 

XIOOO 

Di6 

(nun) 

H50 

(win) 

^84 

(nun) 

Croin) 

SlOl 

5750 

25.5 

0.79 

0.84 

0.15 

0.29 

0.44 

9.0 

S102 

9800 

36.0 

0.95 

0.74 

0.15 

0.29 

0.44 

8.1 

S103 

10515 

43.2 

0.79 

0.47 

0.15 

0.29 

0.44 

8.8 

S104 

10538 

40.5 

0.96 

0.77 

0.15 

0.29 

0.44 

10.7 

S105 

11710 

51.4 

0.80 

0.36 

0.15 

0.29 

0.44 

11.3 

S106‘ 

11858 

41.2 

0.96 

0.62 

0.15 

0.29 

0.44 

10.55 

S107 

12050 

47.8 

0.97 

0.78 

0.15 

0.29 

0.44 

11.2 


Table 4*5 J Hydraulic and sediment data for ribbons on S2 sand bed 
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Table 4.6 : Hydraulic and sediment data for ribbons on S3 sand bed 
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CHAPTER 5 


ANALYSIS OF RESULTS AND DISCUSSIONS 


5.1 OC flow with suet ion/in jaction 

The analysis of OC flow with suet ion/ in j ect ion developed in 
Chapter 3 is semi -empirical in that it uses eddy viscosity 
formulations and assumed criteria for the apparent origin of 
turbulent velocity distribution, which need direct or indirect 
confirmation by experimental data. In the following sections the 
analytical results have been compared with experimental OC and TBL 
data, and significant aspects of the analytical model have been 
discussed in the context of wall flows, in general, and OC flows, 
in particular. 

5.1.1 Free surface elevation 

As described in Chapter 4 twelve sets of OC experiments 

were conducted on nominally rigid, hydraulically rough beds -- six 

sets with influent seepage (suction), and six sets with effluent 

seepage (injection). Each set consisted of an uniform subcritical 

flow (discharge =0© depth = Hj^) run, “and another run 

with suction/injection for the same bed condition and same Q (=0^) 

and free surface slope (S=Sq) at a predetermined location (X=0) . 

Determining the roughness scale. Kg, from an assumed log-law for 

the uniform flow experiment and holding it constant for the 

spatially varied flow, H(x) was computed from the model for 

‘ 2 

-2600< X< 2600. For the present experiments |v„l<<l and (0fg) <<1 
in most cases, and the mixing length (designated a = 2) and 

parabolic eddy viscosity (designated a = 1) models gave 

numerically close solutions. 
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Figure 5.1 shows plots of observed versus computed cha.nge 
in flow depth, AH at X == 0 due to suction/ in j ection. 
The estimated error in measuring A H is indicated in the figure. 
It is seen from the figure that the analytical results for both 
a=l and a=2 show fair to good agreement with experimental results 
within the limits of experimental error, though the mixing length 
model slightly overpredicts |AH| in the case of suction. 

A comparison of the observed and computed water surface 
profiles (elevations with respect to H^) is shown in Figures 5,2a 
and '5.2b for suction and injection respectively. The plotted 
points are the mean values from a set of measurements with 
slightly varying surface profiles, and the measurement error is 
expected to be small. The computed profiles are almost identical 
for a=l and a=2 ; hence only one set of profiles is shown in the 
figure. The figure shows good agreement with experimental data 
for both suction and injection. 

A comparison of the present model with one-dimensional OC 
analysis is in order. For the flow situation studied here the 
friction slope, Sf, used in OC hydraulics is easily determined. 
Setting dH/dX = 0 in Equations 2.36 and 2.39, 


ff 



^ ^m ^w '^av ' 

Vy > 

0 

5.1a 

So 


1 ' 

^e ^w '^av’ 

V 

> 

0 

5.1b 


Uith an assumed 1/7-th power law for u(y) , C^, ~ 1.015 and ~ 
1.045. C« and C.. which depend on the shape factor of u(y), could 
change slightly depending on the exact distribution of u(y j , but 
the above values are adequate for comparison with the present 
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Figure 5,2 t Comparison o£ observed and aomputed water surface 
profiles t (a) suction ; (b) Injection. 


Table 5.1 » Computed parameters for self-similar smooth bed OC 
> • ‘ flow with suction/in Jection * parabolic eddy 
viscosity model (a=l) . 
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Table 5.2 s Computed parameters for self-similar smooth bed 
OC flow with suction/imjection : mixing length 
model (a=2) . 
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Table 5,3 * Computed parameters for self-similar OC flow 
over hydraulically rough bed with suction/ 
injection : k = 100, a = 1. 
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Table 5.4 t Computed parameters for self-similar OC flow over 

hydraulically rough bed with suction/injection * 

“k = 1000, a = 1. 
so 
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diacharc© is roughly proportional to Equations 5.1 give 

Hq _ ' '^w > ° ^-2* 

Hn I Cl - 1.05 < 0 5.2b 

Equation 5.2a agrees fairly with both experimental and analytical 
results, but Equation 5.2b systematically underpredicts |AH|. 

The comparison between Equations 5.2 and the present theory 
can be more extensively made with reference to Tables 5.1, 5.2, 
5.3 and 5.4 which give the flow parameters computed from the 
present model. Tables 5". 1 and 5.2 give the computed parameters 
for smooth bed flow with a=l and a=2 respectively, while Tables 
5.3 and 5.4 are for hydraulically rough bed with kg^ = 100 and 

ro 

1000 respectively, where kg^ = Kg and with a=l 

Between Tables 5.1 and 5.2 the computed parameters are seen to 

diverge with increase in |Vy| . This difference may be due to 

approximations in the analysis for a=2 wherein terms of order 
2 3 

(0£g) I fg) » etc. where ignored. As evident from the magnitude 

2 

of in Table 5.1, (0 f g) may became significant even at Iv^l = 

0.005 since 1 T^-1 | is not negligible. Hence in the rest of this 
chapter attention will mainly be focussed on the parabolic eddy 
viscosity model. 

Comparison of Equations 5.2 with T or Sf/S^ given in the 
above tables shows that the momentum equation generally gives 
comparable results for both suction/injection, whereas the energy 
equation yields lower values of (S£/Sq - 1| in both cases. It is 
therefore suggested that the application of energy equation (with 
the Assumption of constant energy per unit mass) to OC flow with 
suet ion/ in j ect ion may be incorrect, and the momentum equation is 
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recommended for both suction and injection when one dimensional 
analysis is used. This agrees with the recommendation by some 
earlier researchers [cf. El-Khashab and Smith, 1976]. It may be 
noted here that the one-dimensional analysis does not give any 
ready estimate of and dS£/dX near X=0. Conventional one— 
dimensional analysis, wherein dS£/dX is assumed to be zero, will 
indicate for both suction and injection that d^H/dX^ is negative 
for subcritical flows and positive for supercritical flows if Cg = 
^m ~ ^ ' whereas the present analysis indicates negative curvature 
for • injection and positive curvature for suction. The 
experimental results support the present model. 

5.1.2 Shear stress and mean velocity distributions 

The present analysis indicates that the velocity profile of 
u(y) is fuller for suction and more curved for injection. 
However, in the major part of the flow uCy) ie higher for suction 
and lower for injection than the corresponding log-law velocity 
for R* . These results qualitatively agree with the sample plots 
of uCy) with suction and injection given by Brink and Oldenziel 
[1974]. However, as mentioned in Chapter 2, their data are not 
adequate for comparison with the present results. 

In order to elicit further empirical confirmation, computed 
values of U(y) have been normalized with U* and plotted against y, 
where y = yU* /V , in Figures 5.3a and 5.3b for two illustrative 
values of the parameters R*, R* = 300 and 10000 respectively, with 
hydraulically smooth bed conditions. The figures show that u(y) 
systematically deviates from the log-law with injection and 
suction. The deviation is positive for injection and negative 
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for suction, and in both cases increases with increase in *y . 
These trends are qualitatively supported by the data of Maclean 
[1991] for flume experiments with abstraction through a perforated 
bottom. However, direct comparison with Maclean’s results are not 
possible for the following reasons : 

(i) Maclean’s data are for v^/Ug^y > 0.04, ie, for relatively 
high values of IVyl, for which the present analytical 
results do not hold. 

(ii) The length of the porous zone in Maclean’s experiments was 
probably inadequate for full development of the inner 
boundary layer. In fact S and Sq were probably very 
different at the measuring station [cf. Maclean and 
Uilletts, 1986]. 

(iii) U* was determined by Maclean from an assumed relation 
between the spatially varied and uniform flow velocities at 
a reference height. 

In the absence of adequate OC data, Stevenson’s (or 
Dorrance and Dore’s) velocity distributions for TBL flow with 
suction/injection have been compared with the present results in 
Figures 5.3a and 5.3b. Though the accuracy of Stevenson’s 
equation is not accepted by all researchers [eg, Kays, 1972; 
Schetz and Nerney, 1977], it has been advocated by many others 
[eg, Schlichting] as it shows reasonable agreements with empirical 
data. As evident from Figures 5.3 u(y) for TBL flow deviates 
significantly from the log-law in the same manner as for OC flow. 
However the divergence from log— law is much more in the case of 
TBL flow. It may be pointed out here that Figures 5.3a and 5.3b 
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Figure 5»3 i Mean velocity distribution scaled with wall 
variables and comparison with TBL data 
(b) R* = 10,000* 







This result is also true for other values of R* and , indicatinc 
that the assumption of identical bed shear stress for the laminar 
and turbulent velocity distributions used in the present analysis 
is almost equivalent to assuming equal stresses at the interface. 

It may be mentioned that the properties (deviations from 
uniform flow, trends, etc.) of the flow parameters noted above are 
true also for rough turbulent flows save that : 

(i) the deviation from the uniform flow case decreases with 
increase in Kg/H^ for constant R* and v^; and 

(ii) the apparent origin of the velocity, y^, and y^ remain 

invariant for fully rough bed flows if KgU*/i^ > 70, 

unless Kg/H^ changes. 

On the basis of the comparisons in this and the previous 
section, it may be concluded that the premises and results of the 
present analysis are reasonably correct. 

5.1.3 Physical significance for turbulent vail flows 

The results cited and discussed in the previous section 
indicate that the effect of suction/injection on different types 
of wall-bounded flows should be, at least qualitatively, similar. 
Measurements reported in literature indicate that U(y) tends 
towards the laminar velocity profile in TBL, duct and OC flows 
when subjected to injection. The present work confirms this 
result and also indicates that the location of the sublayer- 
turbulent interface, y^, shifts away from the wall in smooth-bed 
flows subjected to injection (see Tables 5.1 and 5.2). Conversely 
U(y) tends to more uniformity and the turbulence region extends 
closer to the wall with suction. Moreover, T (y) increases/ 
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decreases with suction/injection. On the face of it, these 
changes seem to suggest that suction/injection enhances /diminishes 
the turbulent nature of flow, which is contradicted by the 
observed reduction in flow turbulence reported in literature 
[refer eg. Elena, 1984; Fulachier et al , 1987]. The apparent 

conflict is resolved by observing that while auction/injection 
reduces/increases and/or it may vastly increase/ 

decrease the correlation between U’ and V’ . This will cause 
increase in the turbulent shear stress with simultaneous decrease 
in turbulent normal stresses and turbulence generation with 
suction, and the converse for injection. Thus, suction can still 
be considered to promote the turbulent nature of the flow 
(increased ratio of turbulent to "laminar” shear stresses) at the 
cost of less turbulence or with apparent ”laminarizat ion” of the 
flow. This suggests that suction is likely to reduce the strength 
(or vorticity) and/or number of turbulence eddies while enhancing 
their directional preference and reducing their randomness. The 
opposite effects are implied for injection. 

The present analysis makes use of eddy viscosity 
formulations based on Prandtl’s approach. Theoretically, there is 
little to choose between any of these formulae as they are based 
on inadequate postulates or phenomenological hypotheses. It may 
also be mentioned that Prandtl’s mixing length formula as well as 
the two eddy viscosity formulations used in this work give closely 
similar velocity distributions for the inner region of canonical 
TBLs that are well within the scatter of experimental data. 
Physically, Prandtl’s formulation is based on the proportionality 
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of u’ and v' with ydu/dy. The parabolic V formula can be 
interpreted differently by aesuming v' to be independent of du/dy. 
Empirical data [see eg. Hinze, 1975; Sabot et al , 1977] indicates 
that in many canonical wall flows 

V’rms - Cl U*Nr 1 -C 2 y 5.3 

where is of order 1 and C 2 is 1. C 2 is considerably less 
than 1 in duct flows, which may be due to turbulent fluctuations 
crossing over from one half to the other half of the duct. 
Nonetheless, the component of V’ that correlates with U' may be 
expressed on the basis of Equation 5.3 as 

V’ = U* <TTy 5.4 

Since V’ is mainly caused by ejection of wall fluid elements 
during turbulence bursts it is reasonable to consider V’ as 
independent of U and instead scaling on U*. On the other hand U’ 
will be coupled to Y dU/dY with )c as a proportionality constant 
as in Prandtl’s formulation. Thus the parabolic eddy viscosity 
formulation is readily obtained. 

The above formulation suggests that any small change 
(perturbation) in wall variable is likely to be transmitted across 
the boundary layer by turbulence eddies with a mean vertical 
velocity approximating Equation 5.4. The time taken for the 
effect of the variable to permeate the flow depth will be 2H/U». 
The mean distance, , traversed by the flow during this time is 



5.5 
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or, = 2Ug^y/u*, where X^j nay be considered as the transition 

zone or the developing length of the inner boundary layer, U* 
being the modified shear velocity after introduction of the wall 
disturbance/change. Similarly an inward growing layer can be 
visualized to develop from the boundary layer edge due to change 
in outer/f reestream variable, and derived as in Equation 5.5 
with the difference that U* is taken as the upstream shear 
velocity since the outer disturbance propagates inward through 
interaction with wall or inner region eddies. As an illustration. 
Table 5.5 gives the computed x^j for smooth bed OC flows with 
suction/injection [refer Table 5.2]. It is evident from the table 
that x^ increases/decreases with suction/in j ection and with 
increase/decrease in R«. Similar results can be computed for 
hydraulically rough flows with suction/injection from Tables 5.3 
and 5.4. The computed values of x^ are comparable to those 
reported in literature for TBL flow [Simpson, 1971]. 


Table 

5.5 

: Computed developing 
suction/in j ection 

1 engths , 

Xj/H, for OC 

flows with 

R* 


csi 

O 

o 

1 

-0.01 

0 

0.01 

0.02 

300 


33.0 

33.65 

34.5 

36.0 

37.9 

1000 


37.7 

38.8 

40.5 

43.2 

47.5 

10000 


46.8 

48.6 

52.1 

58.5 

— 
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As in the case of suction/in jection, can be computed for 
change in other vail variables if and u* are known. For 

example, with a step jump/reduction in wall roughness is 

expected to decrease/increase and U* is expected to increase/ 
decrease. Hence the transition from a smooth to rough wall 
boundary layer is likely to be much faster than the opposite 
transition, which agrees with experimental obervations [cf. Smits 
and Uood, 1985]. Equation 5.4 can also give an estimate of the 
dimensions of the inner layer as it developes if U(y) is known. 
As an approximation, if U(y) is nearly uniform, U(y) ~ ^^av' then 


u* X u* X 

hi ~ — - Cl ) 


5.6 


2u. 


4u. 


‘ av ’ “ av 

which shows that the thickness of the internal layer, hi, varies 
quadratically with x. 

Taking the above arguments further it may be reasoned that 
when a wall input is small such that 

(i) ll(y) is nearly self-similar when scaled with either U* or an 
outer flow velocity scale, and 

(ii) the fractional change in U(y) is comparable to the 

fractional change in U* over the length x^j , 

then the flow is nearly self-preserving not only in terms of its 
mean flow but also with respect to coherent structures and the 
turbulence field. This follows from the unique relation between 


the mean and turbulent flow fields underlying parabolic V 


t • 


Thus, in the case of suet ion/ in j ection applied at low rates to OC 
flow, the present work gives the relations : 



lOZ 



Since the mean flow is nearly self-similar for I'V^I << 1. 

2 

U (Xjj.y) U* Cx<iD 2v „ Uav ^ ^ 

U (o ,y) U* (o) u*(l-2v„u^^) 

If, and only if, the right side of Equation 5.7 is numerically 
significant, then it implies different scales of evolution of U’ 
and V’ ie, the coherent flow structures may change significantly 
in the flow direction. Though the present arguments are 

speculative in nature, and it does not appear possible to 
demonstrate from the turbulence energy equation that the flow will 
remain nearly self-preserving in all turbulence moments and at all 
scales of turbulence [in the sense discussed by George, 1989] if 
the right side of Equation 5.7 is numerically much less than 1, 
this result follows from the phenomenological premises of 


parabolic formulation. Similar results as Equation 5.7 can be 
derived for gradually varying (quasi-uniform) or slowly developing 
wall turbulent flows (under the influence of small variations in 
wall inputs) if a nearly self-simiar mean velocity solution can be 
found with an appropriate eddy viscosity model. The wall inputs 


include wall heating since theoretical considerations 


(Reynolds analogy) ‘and experimental data [eg. Fulachier and Dumas, 
1975; Fulachier, et al , 1987; Kays, 1972] suggest that the 
mechanism and dif f usivities of heat and momentum transfer are 


comparable . 
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5.1.4. Practica.1 imp 11 cat Iona for alluvial channels and duct 
flows 

Moat unlined canals and natural channel are subjected to 
infiltration and seepafie losses through the channel boundary over 
long stretches and long periods of time. The rate of 
infiltration, depending on the difference in head of channel water 
with ambient ground water table and the soil porosity, may vary 
with time, but are generally small. The pr'esent anaysis shows 
that in such cases u(y) will tend to be more uniform and 
higher than the corresponding log-law velocity. Since, even small 
changes in Vy of the order of 10 ^ can change Sf /S q significantly, 
particularly for kg << 1, the use of a constant friction slope, or 
Manning’s ’n’ , or Chezy’s ’C’, can give erroneous results if the 
seepage rates vary with time. If there is mass efflux through the 
channel boundary — this may occur for instance, due to excess 
bank storage or due to groundwater recharge by prolonged 
precipitation in a catchment or by over-irrigation — the channel 
water level may be significantly higher than predicted. This can 
compound other problems such as the likelihood of a channel to 
overflow its banks during the passage of flood waves. 

Alluvial channels have erodible boundary and carry sediment 
load. Uhile the present analysis shows that infiltration 
increases the time— average bed shear stress, this does not 
necessarily imply that grain mobility and sediment load are 
enhanced. Firstly, the present work gives the mean Xy, but not 
the fluctuating component of Ty, which is probably more 
important in grain erosion and motion. Secondly, grain erosion is 



etrongly affected by hydrodynamic lift force. Eddy viscosity 
models do not give any estimate of the near-wall pressure or lift 
force on surface grains. The change in lift force is probably 
more strongly related to the vertical seepage pressure gradient 
than the OC velocity fluctuations. Thus the combined effect of 
suction on bed grains would be an increase in drag force and 
reduction in lift force which have opposite effects on grain 
erodibility. Going by empirical data [eg. Brink and Oldenziel, 
1974], it seems fairly certain that suction decreases grain 
velocity and bedload in a channel, while injection or effluent 
seepage has the opposite effects. 

As noted in the earlier section the water surface profile is 
convex with effluent seepage and concave with influent seepage. 
Since |Vy| << 1 in most channels, the curvature of the water 

surface would be barely noticeable in field situations. However, 
over long stretches, the curvature effect may change the flow 


depth and 

Uav 

considerably, 

provided S 


So- 

One 

possible 

consequence 

of 

this over a 

long period 

of 

time 

may 

be 

some 

adjustments 

of 

channel geometry such as 

channel 

width 

or 

bed 

slope. It 

may 

be noted that 

many graded 

or regime 

channels 


exhibit exponentially decreasing bed slope over long stretches. 
Several possible causes have been attributed to this phenomenon 
[see eg.. Garde and Raju, 1985]. Another possible reason 
suggested here is long-term channel adjustment of bed slope 
corresponding to the concave water surface produced by 
infiltration from the channel into the bed. Thus if bed slope, S, 


changes as 
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S - exp(- 2Cx) 
then ^bed “ ^bed,o ^ 

which shows that bed elevation, changes at the same rate as 

water surface elevation if C = - a = - (cf. Equation 3.49). 

Hence the flow will be apparently uniform as the flow depth 
remains constant if V^/(Ug^^jj) is constant along X. Such channel 
adjustment leading to nearly self-similar flow may be a natural 
mechanism for channels to adjust to nonuniform flow. The 
suggestion made here is tentative, and obviously needs further 
examination based on field data and a valid variational principle 
of channel adjustment. 

A direct extension of the present analysis to flow in ducts 
is possible if eddy viscosity formulations are assumed to be valid 
for the entire flow region. Replacing S by 8 / 3 x in the OC 
flow equations, the results for flow in rectangular ducts with 
uniform suet ion/ in j ect ion through opposite boundaries will be 
similar to those for OC flow. Some quantitative difference will 
arise due to the flow depth being constant in duct flow. However, 
it may be expected that, in order to carry a certain Qq at a given 
section, the centreline pressure gradient will be significantly 
higher/lower for injection/suction as compared to the case for 
= 0. The ratio of pressure gradients should roughly be about 

(Hq/Hjj)^ in OC flow for |Vy| << 1. Moreover, gP/ax can be 
qxpected to decrease/increase over x at the rate a = Vy/u^y On 
the other hand, the wall shear may increase/decrease slightly with 
suction/injection but will decrease/increase along x with 
suction/injection. Qualitatively similar results can be expected 



for axiaymmetric duct flows. It may be pointed out here that 
Eckert and Rodi’s measurement [refer Narasimhan and Sreenivasan, 
1979] showed clearly an increasing trend for C£ with x in 
turbulent pipe flow with injection. This partly confirms the 
above observations. However, a surer confirmation could be 
forthcoming from direct measurements of the centreline pressure 
gradient which, as of now, does not appear to have been reported 
in literature. Lacking such confirmation the observations made 
herein on duct flows should be considered tentative, or at least 
approximate, since the application of parabolic eddy viscosity 
model does not seem to be accurate even for canonical (ie, 
uniform) duct flows. 

5.2 Hydraulics of sand ribbons 
5.2.1 Scaling of lineations 

The experiments described in the previous chapter showed 
that SI sand beds comprising fine to medium sands developed bed 
ripples at a little above threshold conditions, while the 
relatively coarser S2 and S3 sand beds produced longitudinal 
bedforms or ribbons for a wide range of flow conditions. For SI 
sands at just above incipient state fine striae were observed in 
some cases. On scaling the lateral spacing, L^, of these 
striations with wall variables it was found that they are of 
dimensions on the order of sublayer streaks. As indicated in 
Table 5.6, d^Q < 6.0, and the bed may be considered hydraulically 
smooth. The bed roughness Reynolds number d 5 o = D 50 was 

computed after correcting the bed shear velocity for sidewall 

cs> 

effect [Vanoni, 1975]. The measured values of l^, varying between 
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~ 130 to 200, are higher than sublayer streaks spacings, but this 
could be due the striations not developing at all streak 
locations, or due to all striations not having been clearly 
observed during the measurements, or due to striations being 
related to buffer layer streaks which have higher spacing as 
compared to sublayer streaks [Smith and Metzler, 1977]. On the 
whole, the striations appear to conform to the low-speed streak 
model of parting lineations and bedload streaks. 


Table 5.6 : Lateral spacing of bed striations scaled with vail 

variables 


Run 

SlOl 

S102 

S103 

S104 

S105 

S106 

S107 

dsO 

5.20 

4.77 

4.98 

5.09 

4.65 

4.60 

5.52 

<*84 

7.88 

7.24 

7.55 

7.72 

7.06 

6.98 

8. 38 

ro 

161 

133 

155 

188 

181 

167 

206 


The sand 

ribbons 

on S2 

and S3 sand 

beds 

were more 

well- 


marked and stable structures. Scaling with wall variables did not 
reveal any specific property of the ribbons. When normalized with 
flow depth, however, was found to be fairly constant with an 
average 1^ = ~ 2.03+0.1. This value differs from some 
values reported in literature suggesting 1^ ~ 3 to 4 , particularly 
for ribbons observed in natural settings. This difference may be 
due to the confining effects of sidewalls; but, as the aspect 
ratio of flow was high (about 8 to 16) it is unlikely that 
sidewalls may have restricted the lateral spacing of ribbons. 
Overall, the consistently uniform value of scaled on flow depth 
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confiraia the secondary flow model for sand ribbons. 

The vertical distension, Ly, of sand ribbons does not appear 
to scale with either wall or outer variables. Ly<<L 2 . and in 
general Ly is about 2 to 3 times the median grain diameter of sand 
ribbons or median bedload grain size. Since bedload motion was 
observed during the experiments to occur mainly over and along the 
edge of ribbons, Ly probably represents the thickness of the upper 
two or three layers of moving bedload. Ef fluent /inf luent seepage, 
which had no perceptible effect on , appeared to 

increase/decrease Ly slightly. 

5.2.2 Bed shear stress distribution 

The analysis of secondary flow in Chapter 3 assumed a 
specific type of roughness variation, and corresponding 
sinusoidally varying components of U, U* and V, 

In y^ = In y^^, [1+e cos(kz)] 
u = u^Cy) + ujCy) cosfkz) 

u* = 1 + u*j^ cos(kz) 

and v(y) = aQy®(l-y)*^cos(kz) 

The solution set for uj^(y), u*j^ and a^ depends on the assumed 

values of m, n, and where relates Uj^(l) to uj^fy^), uj^Cl) ~ 

C^UiCyo)- Numerical solution of Equation 3.74 showed that while 
Sq and u^Cy) are significantly affected by Cj^, ra and n, u*]^ 

depends mostly on £ and k, and to a lesser extent on In Vom* m 

and n. Since the analysis is only approximate, and since the aim 
here is to obtain a rough estimate of “the bed stress distribution, 
variations in . m and n may not significantly affect the use of 
the model results in the present study. 
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In order to calibrate the model the data of Muller and 
Studerua [1979] and Studerus [1982] were used. In their 

experiments with simulated bed ribbons, 60 mm wide roughness 
strips were placed at 160 mm centre to centre having 100 ram wide 
smooth zones between the strips. The flow depth being 80 mm, Zq= 2 
or k=7T in their experiments. Replacing their bed roughness by an 
equivalent sinusoidally varying roughness such that the cross- 
sectionally averaged mean velocities match, their velocity and 
shear stress measurements can be compared with the analysis for 
different assumed , m and n. As an illustration some 

comparisons are shown for Cj^ = 0, m = 2 and n = 4 in Figures 5.6 
through 5.8. Figure 5.6 shows fairly good matching between the 
observed and computed secondary flow velocities. The isovels for 
mean streamwise velocity shown in Figure 5.7 also indicate 
reasonable agreement for y<0.4 but diverge from measured values in 
the upper half of the flow region. This can be corrected by 
assigning a negative value to C;]^ . Figure 5.8 shows a plot of 
measured and computed shear stress over the rough and smooth 
zones. Though the measured and computed stresses have similar 
trends of variation over depth the ccomputed stresses are much 
less above the smooth zone. This can be rectified partly, but not 
entirely, by adjusting , m and n. 

It may be noted that in the experiments of Muller and 
Studerus, 6 ~ 0.2 is not very small. Hence .ignoring terms of 
order £ in the analysis may have produced significant error for 
this flow situation. Partly for this reason and partly because 
the main aim here is to estimate the bed shear stress variation 
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Figure 5«6 t Comparison of observed and computed 
secondary flow velocities on bed 
with roughness strips. 
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Figure 5*8 t Comparison of observed and computed 
shear stress over smooth and rough 
strips 


for small values of ^ , accurate calibration of the model 
parameters have not been attempted. Since the sand ribbons 
generated spacing of about 2 Hq, and were computed for 

different values of e and In with k=2 and , m and n as 

above. The computed values are shown in Table 5.7. 




Table 5.7 Computed values of u*]_ and a^ for flow over sand 
ribbons with k=2, 0^=0, m=2 and n=4. 


y om 

0.01 

0.005 

0.002 

0.001 

0.0005 

^om 

-4.605 

-5.298 

-6.215 

-6.908 

-7.601 


268 

213 

180 

167 

157 

u*i/ e 

-2.88 

-2.29 

-1.92 

-1.78 

-1.67 


5.2.3 Bedform reftine criteria and sand ribbons 

The experiments with SI, S2 and S3 sands showed that sand 
ribbons formed readily for flow with equilibrium sediment load 
over a non-uniform sand bed comprising medium to coarse sand or 
sand-gravel mixtures, but not with fine-medium sand mixtures. 
Several bedform regime criteria have been proposed in sediment 

transportation literature, eg., Garde and Albertson, Simons and 
Richardson [refer Garde and Raju, 1985; Vanoni , 1975], van Rijn 
[1984]; several other criteria are in use in geophysical and 

sedimentological analysis. Since all these criteria are 
empirical, it will suffice to compare the sand ribbons’ existence 
conditions with only one of these in order to show that the 
bedform regime criteria are not valid for sand ribbons. 

Figure 5.9 shows a plot of the regime data of sand ribbons 

on Simons and Richardson’s graph of stream power ‘^w'^av versus 

fall diameter of ribbon sand, where is the estimated stress 

over ribbons. As evident from the figure the sand ribbons occupy 
a considerable area of the lower (dune) and transition bedform 
regime. Likewise, comparison of the hydraulic regime of sand 


Likewise 


11 ^ 


ribbons with other bedform regime criteria indicate that with non- 
uniform sand beds containing sufficient amounts of coarse or very 
coarse sand or gravel the development of ripple or dune bedforms 
are considerably inhibited. Physically, this may be attributed to 
possible mutual cancellation of three-dimensional or streamwise 
disturbances associated with suspended motion (of fine and medium 
sands), saltation (of coarse sand) and bedload motion (of very 
coarse sand and gravel). Since the present experiments where 
conducted at subcritical flow conditions, the results complement 
Chiew’s [1991] finding of no antidune formation on non-uniform 
sand bed at near-critical or supercritical flow conditions. 

5.2.4 Sediment transport with sand ribbons 

Since sand ribbons are flow-aligned and are of low relief, 
they do not offer any significant form resistance. Hence the 
equilibrium sediment transport rate, , over a bed of sand 
ribbons may be expected to differ from that over a dune or ripple 
bed. Though the suggested bedload transport equations for sand 
mixtures are different from those of uniform sand, the Einstein- 
Brown relation [refer Vanoni , 1975] for uniform bedload size has 
been used here as an yardstick for comparison. Accordingly, 




5.8 


5.9 


5.10a 



stream power (Ib/ft/s) 


AA. * A Transitio 






115 


and 


36 P “y ^ 


5.10b 


have been plotted in Fi£ure 5.10. As evident from the figure some 
of the data of S3 sand indicate lover bedload than predicted from 
the Einstien-Brown relation for > 0.09, 


40 


0 . 


5.11 


There ia also slight increase in bedload with injection and 
decrease with suction (not shown in figure). This difference can 
be rectified by including a lift force component in Equation 5.9 
to account for change in sediment mobility with suction/ injection 
[cf. Brink and Oldenziel, 1974; Maclean and Willetts, 1986]. 


5.3 Concluding remarks 

The present work shows that eddy viscosity or mixing length 
models can be successfully applied to OC flows where the secondary 
velocity components, V and U, are numerically much less than the 
streamwise velocity, U. In such cases UCx,y,z) can be split into 
separate components 


uCx,y,z) ~ UqCv) + f(x) ui(y) + gCz) U 2 Cy) 5.12 

where u^Cy) is the canonical (log-law) velocity, and fCx), gCz), 
U]^(y) and U 2 Cy) are determined from the equations of motion with 
suitable assumptions on the form of f(x) and g(z) from boundary 
conditions. Evidently, the approach ia likely to fail when the 
secondary velocities or velocity gradients are significant in 



comparison to the streamwise velocity or velocity gradient, eg., 
flow near corners or sidewalls in a rectangular channel, or 
spatially varied OC flow with high rates of inflow or outflow. 
Though the analysis of TBL flow, spatially varied pipe flow and 
other slowly developing wall flows is beyond the scope of this 
work, it is suggested that eddy viscosity modelling together with 
the assumption of Equation5,12 may be found analytically useful in 
many such flow situations. 



CHAPTER 6 


SUMMARY AND CONCLUSIONS 


6 . 1 Summary 

Isotropic eddy viscosity modelling of Reynolds stresses 
has been applied to open channel flow 

(i) subjected to uniformly distributed suction/injection, and 

(ii) over sand ribbons of lateral grain sorting. 

In the former case the equations of motion are reduced to 
two coupled ordinary differential equations by assuming the mean 
streamwise velocity and the outer and inner velocity scales 
(energy gradient velocity and bed shear velocity) to be piecewise 
linear in X. For low rates of suction/injection and water surface 
nearly parallel to the bed, the flow was found to be nearly self- 
similar, and the velocity and shear stress distributions were 
obtained. The analysis indicates higher/lower mean velocity 
(hence lower/higher flow depth) and positive/negative water 
surface curvature with suction/injection as compared to an 
equivalent uniform flow. Flume experiments conducted on rigid, 
hydraulically rough bed showed good agreement with the model. The 
results also compare with boundary layer data while indicating 
quantitative differences. 

Experimental study of sand ribbons showed stable flow- 
aligned sand ribbons at the exclusion of ripplefi and dunes in 
subcritical flow with equilibrium sediment load on a bed of sand- 
gravel or fine-coarse sand mixtures. The spanwise spacing of 
ribbons was about 2 times the flow depth and their relief was 

about twice the median grain diameter. The analytical model 
indicated secondary flow cells and spanwise shear stress variation 
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at the rate of 2 to 3 times the rate of log-law velocity shift due 

to bed roughness variation. 

6.2 Conclusions 

The following conclusions are drawn in the present study : 

2 

Ci) U* /CagHoS^), or Sf/S^ where S£ is friction slope, 
increases/decreases with suction/injection in agreement with 
one-dimensional OC analysis and TBL data. 

(ii) CC and 3 are comparable, but deviate as |Vyl increases. 

(iii) Dimensionless velocity u(y) increases and has a fuller 
profile for suction, and conversely for injection. 

(iv) T(y) sharply increases/decreases near the bed with 
suction/ inj ect ion. 

(v) The above changes magnify with increase in Hq 
and with decrease in relative roughness, Kg/H^. 

Cvi) Water surface has concave/convex curvature for 
suet ion/ in j ect ion . Flow depth varies as 

H ~ Hq [1 - a SqX^ {1 + ' Ca / 3 Dx) ] C3) 

(vii) The analysis suggests that in duct flow dP /dX 
decreases/increases with suction/injection for same 
discharge . 

(vili) From a physical interpretation of parabolic eddy viscosity 
distribution it is reasoned that the flow is nearly self- 
preserving in all respects -- mean flow, turbulence field 
and coherent structures — if lv^^U/U*| << 1, and an 

estimate of inner boundary layer length with change in wall 


inputs is suggested. 
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(ix) It ia aufigeated that affluent aeapafte incraaaaa flow dapth 
and anhancaa chancaa of flood wavaa ovarflovlna atreaiw 
banka, while infiltration may be the cauae of exponentially 
decreasing bed slope in firaded streams. 

(x) Sand ribbons, that are hydrodynamicaily similar to flat bed, 

stable phenomenon associated with coarse sand-fine 
sand or sand-gravel mixtures. 

(xi) Sand ribbons and secondary flow cells develop by mutual 
reinforcement, secondary flow causes bedload deflection 
towards the ribbons, while the ribbons themselves produce 
lateral shear stress variation and secondary flow by 
creating transverse velocity gradients, and hence lateral 
shear stress gradients. 

(xii) Sand ribbons do not grow into well-developed ridge-groove 
structures (except under degrading bed conditions) probably 
because the lateral stresses are not very high, and because 
developed ridges will have the opposite effects as grain- 
sorting, thereby weakening the secondary flow and shear 
stress variation. 

(xiii) The breakdown of ribbons to primary ripples or dunes in 
sand mixtures may be inhibited due to insufficiency of fine 
8and( suspended and saltating sediment grains ). 

(xiv) The existence of sand ribbons in subcritical flow defies the 
existence criteria of ripples, dunes, etc., for uniform 
sand, but complements Chiew’s finding of no antidune bed for 
sand mixjtures in supercritical flow. 
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(xv) From the point of view of paleocurrent analysis lineations 
with lateral grain sorting should indicate past flow with 
orin bedload, whereas uniform bedding surface texture 
may be the imprint of sublayer flow structure. 

6.3 Recommendations 

(i) The scope of the present analytical solution for 

auction/ln j ection can be improved by considering arbitrary 
water surface slopes or higher suction/ in j ection velocities. 
This will necessitate numerical solution of coupled 

equations derived from the equations of motion (boundary 
layer equations). 

(ii) The effect of suction/injection on grain mobility and 

sediment transport rates has been studied by several 
researchers empirically. These studies may be combined with 
the present analytical results. 

(iii) An extension of the present analysis to duct flows and TBLs 
may be studied theoretically and experimentally by adopting 
an appropriate eddy viscosity formula. 

(iv) The sand ribbons studied here were generated in the 
laboratory. Field studies may be a useful addition to the 
knowledge of these structures. 

(v) The study of the breakdown of sand ribbon to ripples or 
dunes (associated with flow disturbances propagating 
downstream) may be a stimulating problem that can throw 
light on the mechanism of generation of ripples and dunes. 
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